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Abstract 

We use heat kernels or eigenfunctions of the Laplacian to construct local coordinates on large 
classes of Euclidean domains and Riemannian manifolds (not necessarily smooth, e.g. with C a 
metric). These coordinates are bi-Lipschitz on embedded balls of the domain or manifold, 
with distortion constants that depend only on natural geometric properties of the domain or 
manifold. The proof of these results relies on estimates, from above and below, for the heat 
kernel and its gradient, as well as for the eigenfunctions of the Laplacian and their gradient. 
These estimates hold in the non-smooth category, and are stable with respect to perturbations 
within this category. Finally, these coordinate systems are intrinsic and efficiently computable, 
and are of value in applications. 
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1 Introduction 

The concept of a coordinate chart for a manifold is quite old, but it has only recently become a 
subject of intensive study for data sets. In this paper we will state and prove a new theorem for 
coordinate charts on Riemannian manifolds. This result is meant to explain the empirically observed 
robustness of certain coordinate charts for data sets, Before stating our results, we explain in more 
detail the setting, first for manifolds, and then for data sets. 

Let M. be a Riemannian manifold. A coordinate chart (more precisely, a restriction of one) can 
be viewed as a mapping from a metric ball B C M. into R d , where d is the topological dimension of 
M. This mapping has the form 

F(x) = (f 1 (x),f 2 (x),...J d (x)). 

It is natural to ask for F to have low distortion. Let F(B) = B C M. d . By assumption F is a one to 
one mapping from B to B. The Lipschitz norm of F is defined as 

_ \\F(x)-F(y)\\ 

\\r\\Lip- SUp r 

x,»gB dM\x,y) 

where dj^i'i ') is the metric on M. and || • || is the usual Euclidean metric on M d . Similarly, one sets 

\\F \ L i P = sup — — 

*, y jB \\F{x) - F(y)\\ 

Then the distortion of F on B is defined to be 

Distortion^, B) := \\F\\ Llp x WF^W^p. (1.0.1) 

It is worth recalling at this point a prime example of a coordinate chart, namely the coordinate 
chart on a simply connected planar domain T> given by a Riemann mapping F from T> to the unit disc 
D. Let z € T> and define r = dist(z , <9P). If we choose our Riemann map F to satisfy F(z ) = 0, 
then the distortion theorems of classical complex analysis (see e.g. [30] page 21) state that F maps 
the disc B(zq,^) onto "almost" the unit disc, with low distortion: 

B(0, k- 1 ) C F(B(z , T -)) C B(0, 1 - kT 1 ) , (1.0.2) 

Distortion^, B(z , ^)) < « . (1.0.3) 

In other words, on B(zq 1 |), i^ 1 is a perturbation (in the proper sense) of the linear map given by 
F'(z )(z-z ), and |F'(z )| ~ ± 
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In this paper we will look for an analogue of (|1.0.2p and (|1.0.3[) above, but in the setting of 
Riemannian manifolds. We will show that on Riemannian manifolds of finite volume there is a 
locally defined F that has these properties, and that this choice of F will come from globally defined 
Laplacian eigenfunctions. On a metric embedded ball B C M. we will choose global Laplacian 
eigenfunctions ip^,^^, ■■■,<Pi d and constants 71, 72, 7<z < k (for a universal constant n) and define 

* := (iWh > 72^2 > ■••) Jdfid) • (1.0.4) 

This choice of $, depending heavily on zq and r, is globally defined, and on B(zq, n~ l r) enjoys the 
same properties as the Riemann map does in (j 1 .0. 2() and (|1.0.3|> . In other words, <1> maps B(zo, « _1 r) 
to, roughly, a ball of unit size, with low distortion. Here we should point out the 1994 paper of 
Berard et al.[5] where a weighted infinite sequence of eigenfunctions is shown to provide a global 
coordinate system (points in the manifold are mapped to 1%). To our knowledge this was the first 
result of this type in Riemannian geometry. Our results can be viewed as a strengthening of their 
work, and have as a consequence the statement that for a compact manifold without boundary, a 
good global coordinate system is given by the eigenfunctions tfj with eigenvalues Xj < kR^. Here 
i?inj is the inradius of M, i.e. the largest r > such that for all x € Ai, B(x,r) is an embedded 
ball. 

The impetus for this paper and its results comes from certain recent results in the analysis of 
data sets. A recurrent idea is to approximate a data set, or a portion of it, lying in high dimensional 
space, by a manifold of low dimension, and find a parametrization of such data set or manifold. This 
process sometimes goes under the name of manifold learning, or linear or nonlinear dimensionality 
reduction. This type of work has been in part motivated by spectral graph theory [8] and spectral 
geometry [TJ [24J [18] (and references therein). Let {xj}^ be a collection of data points in a metric 
space X. It is frequently quite difficult to extract any information from the data as it is presented. 
One solution is to embed the points in W 1 for n perhaps quite large, and then use linear methods 
(e.g. those using singular value decomposition) to obtain a dimensional reduction of the data set. In 
certain situations however linear methods are insufficient. For this reason, there has recently been 
great interest in nonlinear methods 0- Unfortunately such techniques seldomly come with guarantees 
on their capabilities of indeed finding local parametrization (but see, for example, [HI [20l [54] ) , or 
on quantitative statements on the quality of such paramctrizations. 

One of these methods, diffusion geometry, operates by first defining a kernel K{xj, Xk) on the data 
set, and then altering this slightly to obtain a self-adjoint matrix (rrij^) that roughly corresponds to 
the generator of a diffusion process. The eigenvectors of the matrix, should be seen as corresponding 
to Laplacian eigenfunctions on a manifold. One (judiciously) selects a collection Vi 1 , Vi 2 , Vj„ . of 
eigenvectors and maps 

Xk -» (uii , v i2 , v im ) G M m (1.0.5) 

Careful choices of collections of eigenvectors have been empirically observed to give excellent repre- 
sentations of the data in a very low dimensional Euclidean space. What has been unclear is why this 
method should prove so successful. Our results show that in the case of Riemannian manifolds, one 
can prove that this philosophy is not just correct, but also robust. It is to be said that researchers so 
far have restricted their attention to the case when the lowest frequency eigenfunctions are selected, 

i.e. ii = i,» 2 = 2,...,i m = m nam no [gin]. 

Given these results, it is plausible to guess that an analogous result should hold for a local piece 
of a data set if that piece has in some sense a "local dimension" approximately d. There are certain 
difficulties with this philosophy. The first is that graph eigenfunctions are global objects and any 

1 Examples of such disparate applications include document analysis [15] . face recognition '21 i . clustering [371 IT], 
machine learning 4 38, 52 34, 33 53 , nonlinear image denoising and segmentation 45, 52_, processing of articulated 
images [20] , cataloguing of galaxies [21] , pattern analysis of brain potentials 32 and EEG data 44 , and the study of 
brain tumors [6]. A variety of algorithms for manifold learning have been proposed 41, 3 4 31, 13 9 10 28 59 20 
I55ll59ll56ll54ll43ll42] . 
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definition of "local dimension" may change from point to point in the data set. A second difficulty 
is that our results for manifolds depend on classical estimates for eigenfunctions. This smoothness 
may be lacking in graph eigenfunctions. 

It turns out that another of our manifold results does not suffer from these serious problems 
when working on a data set. We introduce simple "heat coordinate" systems on manifolds. Roughly 
speaking (and in the language of the previous paragraph) these are d choices of manifold heat kernels 
that form a robust coordinate system on B(zq, n~ l r). We call this method "heat triangulation" in 
analogy with triangulation as practiced in surveying, cartography, navigation, and modern GPS. 
Indeed our method is a simple translation of these classical triangulation methods, and has a closed 
formula on M. d , which we note has infinite volume! (Our result on heat kernels makes no assumptions 
on the volume of the manifold.) For data sets, heat triangulation is a much more stable object than 
eigenfunction coordinates because: 

• Heat kernels are local objects (see e.g. Proposition 13. 3.2]) 

• If a manifold M is approximated by discrete sets X, the corresponding graph heat kernels 
converge rather nicely to the manifold heat kernel. This is studied for example in [3T | IT3 l \14 \ 13] . 

• One has good statistical control on smoothness of the heat kernel, simply because one can 
easily examine it and because one can use the Hilbert space {/ <E L 2 : V/ £ L 2 }. 

• Our results that use eigenfunctions rely in a crucial manner on Weyl's Lemma, whereas heat 
kernel estimates do not. 

In a future paper we will return to applications of this method to data sets. 
The philosophy used in this paper is as follows. 

Step 1. Find suitable points yj, 1 < j < d and a time t so that the mapping given by heat kernels 
(x — > Kt(x,yi), K t (x, yd)) is a good local coordinate system on B(z, K~ x r). (This is heat 
triangulation.) 

Step 2. Use Weyl's Lemma to find suitable eigenfunctions <fi j so that (with Kj(x) = K t (x,yj)) one 
has large gradient. 

Each point y <= M. gives rise to a heat kernel Kt(x,y). One may think of Step 1 as sampling this 
family of heat kernels K t (x,y) at d different choices y\ 1 ■■■,yd- Indeed, with high probability, ran- 
domly chosen points from the appropriate annulus will be suitable. Step 2 corresponds to sampling 
the vector {ipj(x)e Xjt }j d times, once for each point y\, ...,yd- This last sampling, where we choose 
an index j, cannot be performed randomly! (See example in section f5.ip . 

At this point we would like to note an advantage that local parametrization by eigenfunctions has 
over heat kernel triangulation (which we do not discuss in this paper). Consider the planar domain 
[0,3e] x [0,3]. Then, using only two Neumann eigenfunctions, one gets a good parametrization of 
the rectangle [e, 2e] x [1,2]. On the other hand, in order to get parametrization of similar distortion 
using heat kernel triangulation, on needs to use ~ - different heat kernels. 

To see where our philosophy comes from, we return for a moment to the setting of a simply 
connected planar domain T> of area = 1. Let z £ T> and r be as in the discussion before equation 
(|1.0.2[) . With the choice of Riemann mapping F, with F(zq) — we have the classical formula 
known to Riemann: 



Here G(-,z ) is Green's function for the domain T>, with pole at z , and G* is the multivalued 
conjugate of G. Thus, all information about F on B(zq, §) is encoded in G(z, Zq). Recall that 



F(z) = exp{ - G{z, zo) - iG* (z, z ) } . 



(1.0.6) 



X 




(1.0.7) 



n 
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where K is the (Dirichlet) heat kernel for T>. Thus the behavior of F can be read off the information 
on K(z, zq, t). Now write, 

oo 

K(z,z ,t) =^^(z)^(z )e A ^ (1.0.8) 
j'=i 

where {<fij} is the collection of Dirichlet eigenfunctions (normalized to have L 2 norm = 1) and 
Aifj = Xjifj. Notice that 

\F>(z)\ = \VG(z,z )\e- G ( z - z «K (1.0.9) 

Since ~ 1 on B(zo, |) it is reasonable to guess from the above identities that there are 

eigenfunctions ipj such that 

|V^|>i (1.0.10) 

on B(z , K~ 1 r), for some n > 1, independent of T>. (More precisely, a short calculation with Weyl's 
estimates makes this reasonable.) This simple reasoning turns out to be correct and the main idea 
of this paper. The proof does not depend on any properties of holomorphic functions, but runs with 
equal ease in any dimension. This is because it only requires estimates on the heat kernel, Laplacian 
eigenfunctions and their derivatives, all of which are real variable objects. 

The paper is organized in a top-bottom fashion, as follows. In Section [2] we state the main 
results, in Section [3] we present the main Lemmata, the proofs of the main results, and important 
estimates on the heat kernel and eigenfunctions of the Laplacian, together with their proofs, but 
mostly only in the Euclidean case. For the purpose of completeness we have recorded here proofs 
of several known estimates, over which the experts may wish to skip. In Section [3.51 we present the 
material for generalizing most estimates to the manifold case. Finally we discuss some examples in 
Section [SJ We include a Table of notation at the end of the manuscript, see Section [5J 



2 Results 

2.1 Euclidean domains 

We first present the case of Euclidean domains. While our results in this setting follow from the 
more general results for manifolds discussed in the next section, the case of Euclidean domains is 
of independent interest, and the exposition of the main result as well as the proof in this case is 
simpler in the several technical respects. 

We consider the heat equation in fi, a finite volume domain in M d , with cither Dirichlet or 
Neumann boundary conditions i.e., respectively, 

(A-&)u(z,t) = or J(A-£Mx,f) = 

Here v is the outer normal on <9f2. Independently of the boundary conditions, A denotes the 
Laplacian on Q. In this paper we restrict our attention to domains where the spectrum is discrete 
and the corresponding heat kernel can be written as 

K t (z,w) = K?(z,w) = , £<p j (z)ip J (w)e- x ' t . (2-1.1) 

3=0 

where the {(fj} form an orthonormal basis of eigenfunctions of A, with eigenvalues < Ao < • • • < 
Xj < .... We also require a (non-asymptotic) Weyl-type estimate: there is a constant C coun t such 
that for any T > 

#{j : < Xj < T} < CcountTi |fi| . (2.1.2) 
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In the Dirichlet case C coun t does not depend on f2 (see remark [3.4.3P . For the Dirichlet case the 
only substantial problem is that the eigenfunctions may fail to vanish at the boundary. This in 
turn only occurs if there are boundary points where the Wiener series (for the boundary) converges 
[571 I5U] . For the Neumann case the situation is more complicated [57J [35] . In particular there are 
domains with arbitrary closed continuous Neumann spectrum |27j . We therefore restrict ourselves 
in this paper to domains (and, later, manifolds) where conditions (|2.1.ip and (|2.1.2[) are valid. More 
general boundary conditions can be handled in similar fashion, since our analysis is local and depends 
on the boundary conditions only through the properties above. 

Finally, here and throughout the manuscript, we define f f '■= rm J f- 

B B 

Theorem 2.1.1 (Embedding via Eigenfunctions, for Euclidean domains). Let SI be a finite volume 
domain in M. d , rescaled so that = 1. Let A be the Laplacian in il, with Dirichlet or Neumann 
boundary conditions, and assume that ([2.1. 1[) and (|2.1.2|) hold. Then is a constant n > 1 that 
depends only on d such that the following hold. 

For any z £ fi, let p < dist (z, dtt). Then there exist integers ij., . . . , i<j such that, if we let 



, 1 = 1,. ..,d, 

we have that: 

(a) the map 

$:B(z,k-V) -> K d (2.1.3) 

x ' ^ (■y 1 (p ll (x),...,j d ip ld (x)) (2.1.4) 

satisfies, for any x\,x 2 £ B(z, n~ 1 p), 

— \\xi - s 2 || < ||$(x x ) - $(z 2 )|| < -\\xx - sail ; (2.1.5) 
P P 

(b) the associated eigenvalues satisfy 

n~ 1 p~ 2 < X h , . . . , X id < np~ 2 ; 

(c) the constants "fi satisfy 

{C count) 2 • 

Remark 2.1.2. In item (c) above, it will also be the case that k^ 1 pi < jj. 

Remark 2.1.3. The dependence on C coun t is only needed in the Neumann case because, unlike the 
Dirichlet case, the upper bound in Weyl's Theorem depends on the domain. See Remark \3.4-3\ for a 
more precise statement. 

2.2 Manifolds with C a metric 

The results above can be extended to certain classes of manifolds. In order to formulate a result 
corresponding to Theorem 12.1.11 we must first carefully define the manifold analogue of dist (2, dfl). 
Let M. be a smooth, (/-dimensional compact manifold, possibly with boundary. Suppose we are 
given a metric tensor g on M. which is C a for some a £ (0, 1]. For any zq £ Ai, let (U,F) be a 
coordinate chart such that zq £ U and normalized so that 
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(i) g a (F(z )) =S il . 
Then we assume that 

(ii) for any x E U, and any £, v E R d , 

d d 

CnM\\Z\\l* < E 9 ij (F(x)Mj and J2 <^(F(z))6^ < c max ( 9 )||£|| K<i \\v\\ Rd . (2.2.1) 

We let 

r a (2o) = sup{r > : B r (F(z )) C F(t/)} . (2.2.2) 

Observe that, when g is at least C 2 , can be taken to be less than the inradius, with local coordinate 
chart (U, F) given by the exponential map at z. The chart (U, F) may intersect the boundary with 
no consequence, as all of the work will be done inside B(zo,ru)- We denote by \\g\\ a the maximum 
over all i, j of 

\g^(F(x))-g^(F(y))\ 
J y \F(x)-F(y)\<* 

for x, y in U. The natural volume measure dp on the manifold is given, in any local chart, by i/det g ; 
conditions (|2.2.1[) guarantee in particular that det<? is uniformly bounded below from 0. Let Am be 
the Laplace Beltrami operator on A4. In a local chart, we have 

A M f(x) = --/=|== E d i {\f^99 ij {F{x))dif) (F{x)) , (2.2.3) 

when g is smooth enough (e.g. g E C 1 ). In general one defines the Laplacian through its associated 
quadratic form [17l H"6] . Conditions (|2 . 2 . 1 [) are the usual uniform ellipticity conditions for the 
operator (|2.2.3p . With Dirichlet or Neumann boundary conditions, Am is self-adjoint on L 2 (Ai, p). 
We will assume that the spectrum is discrete, denote by < Ao < • • ■ < \j < its eigenvalues and by 
{(fj} the corresponding orthonormal basis of eigenfunctions, and write equations (|2.1.ip and (|2.1.2|) 
with fl replaced by M.. 

Theorem 2.2.1 (Embedding via Eigenfunctions, for Manifolds). Let (A4,g), z E M. be a d di- 
mensional manifold and (U,F) be a chart as above. Assume \M\ = 1. There is a constant k > 1, 
depending on d, c m i n , c maX ; ||<?IU; ol, such that the following hold. 

Let p < rjj{z). Then there exist integers i\, . . . ,id such that if we let 



li 



j rf | . / = 1 d - 

\B(z,K-ip) 



we have that: 
(a) the map 



$:B(z,k~V) -> ^ d (2.2.4) 
x ' ^ (j 1 ip ll (x),...,j d ip ld (x)) (2.2.5) 



satisfies for any X\,x% E B(z, K 1 p) 



d M (x u x 2 ) < \\$(xi) ~ $(X2)\\ < - d M (x u x 2 ). (2.2.6) 

P P 
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(b) the associated eigenvalues satisfy 

k~ x P~ 2 < Ai lf . . . , \ id < np~ 2 . 

(c) the constants 7/ satisfy 

7i. •••>7d < K(C co „„ t )2 . 

Remark 2.2.2. As in the Euclidean case, in item (c) above, it will also be the case that kT 1 pi < -fj 

Remark 2.2.3. Most of the proof is done on the local chart (U,F) containing z. An inspection of 
the proof shows that we use only the norm \\g\\ a of the g restricted to this chart. 

Remark 2.2.4. When rescaling Theorem \2.2.1[ it is important to note that if f is a Holder function 
with \\f\\c = A and f r (z) = f(r~ 1 z) then \\f r \\c a — Ar a . Since we will have r < 1, f r satisfies a 
better Holder estimate then f, i.e. 

\\f r \\ Ca =Ar a <A=\\f\\ c «. 
We will repeatedly use this observation when discussing manifolds with C a metric. 

Remark 2.2.5. We do not know, in both Theorem \2.1.1\ and Theorem \2.2.1\ whether it is possible to 
choose eigenf unctions such that 71 ~ 72 ~ ... ~ 7<j. If this were so, the map x 1— > (tp^ (x), . . . , tp id (x)) 
would be a low distortion map whose image has diameter > k _1 . 

Remark 2.2.6. As was noted by L. Guibas, when M. has a boundary, in the case of Neumann 
boundary values, one may consider the "doubled" manifold, and may apply our result for a possibly 
larger rjj(z). 

Clearly Theorem 12.1.11 is a particular case of Theorem 12.2.11 but the proof of the former is 
significantly easier in that one can use standard estimates on eigenfunctions of the Laplacian and 
their derivatives. For the sake of presentation we present one proof for both Theorems, but two sets 
of required Lemmata for those estimates which are significantly different in the two cases. 

Remark 2.2.7. The method of the proofs also gives a result independent of the constant C coun t'. 
Let (M, g) and z € M. be as in Theorem \2.2.1\ Let n > 0, and assume that for any x G M. we have 
a chart (U,F) such that rjj{x) > n > (in particular, A4 has no boundary). Then for p < r\ the 
same results as in Theorem \2.2.1\ hold, except the constant k depends only on d, c m i n , c max , ||<?||a; 
a and not on C coun t- This is due to the fact that C coun t becomes universal for values of T > w~ 2 . 

Another, in some sense stronger, result is true. One may replace the d eigenfunctions in Theorem 

12.2. ll bv d heat kernels {K t (z, yi)}i=x 4. In fact such heat kernels arise naturally in the main steps 

of the proofs of Theorem 12.1.11 and Theorem 12.2.11 This leads to an embedding map with even 
stronger guarantees: 

Theorem 2.2.8 (Heat Triangulation Theorem). Let (A4,g), z S M. and (U,F) be as above, with 
the exception we now make no assumptions on the finiteness of the volume of M. and the existence 
of C coun t. Let p < rjj(z). Let pi, ...,Pd be d linearly independent directions. There are constants 
c > and c', k > 1, depending on d, c m i n , c max , p Q ||5||aj ct, and the smallest and largest eigenvalues 
of the Gramian matrix ({pi,Pj))i,j=i,...,d> such that the following holds. Let yi be so that yi — z is in 
the direction pi, with cp < dM{yii z ) < 2cp for each i = 1, . . . ,d and let t = kT 1 p 2 . The map 

x 1 ► (p d K t (x, yi)), p d K t (x, y d )) (2.2.7) 

satisfies, for any Xi,X2 € B(z, n~ 1 p), 

^-—d M {xi,x 2 ) < \\®{xi) - $(a;2)|| < —d M (xi,x 2 ) ■ (2.2.8) 
dp p 
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The reason for the factor p a which we have in p a ||g|| a above is to get scaling invariance. 

This theorem holds for the manifold and Euclidean case alike, and depends only on the heat 
kernel estimates (and its gradient). We again note that for this particular Theorem we require no 
statement about the volume of the manifold, the existence of L 2 Laplacian eigenfunctions, or their 
number. The constants for the Euclidean case, depend only on dimension, and not on the domain. 
The content of this theorem is that one is able to choose the directions y, — z randomly on a sphere, 
and with high probability on gets a low distortion map. This gives rise to a sampling theorem. 

One may replace the (global) heat kernel above with a local heat kernel, i.e. the heat kernel for 
the ball B{z, p) with the metric induced by the manifold and Dirichlet boundary conditions. In fact, 
this is a key idea in the proof of all of the above Theorems. Thus, on the one hand our results are 
local, i.e. independent of the global geometry of the manifold, yet on the other hand they are in 
terms of global eigenfunctions. 

As is clear from the proof, all theorems hold for more general boundary conditions. This is 
especially true for the Heat Triangulation Theorem, which does not even depend on the existence 
of a spectral expansion for the heat kernel. 

Example 2.2.9. It is a simple matter to verify this Theorem for the case where the manifold in 
M. d . For example if d = 2, p = 1, and z = 0, y\ = (—1, 0) and j/2 = (0, —1). Then if K t (x, y) is the 
Euclidean heat kernel, 

x -> (Ki(x,yt), Ki(x,y 2 )) 

is a (nice) biLipschitz map on B[(0, 0), i). (The result for arbitrary radii then follows from a scaling 
argument). This is because on can simply evaluate the heat kernel 



So in Si ((0,0)) 



VKx{x, Vl ) ~ -U-i(l,Q) and Vffi(a;,i&) ~ ^-e~i{0, 1) . 

Z7T Z7T 
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3 The Proof of Theorems [2331 and I2T2TT1 

The proofs in the Euclidean and manifold case are similar. In this section we present the steps of 
the proofs of Theorems 12.1. 11 \2. 2. II we will postpone the technical estimates needed to later sections. 

Because we may change base points, we will use R z (or similarly, R w ) in place of p. We will also 
interchange between B(x,r) and B r (x). 

Remark 3.0.10 (Some remarks about the Manifold case). (a) As mentioned in Remark 1 2. 2.3[ 
we will often restrict to working on a single (fixed!) chart in local coordinates. When we 
discuss moving in a direction p, we mean in the local coordinates. 

(b) We will use Brownian motion arguments (on the manifold). In order to have existence and 
uniqueness one needs smoothness assumptions on the metric (say, C 2 , albeit less would suffice, 
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see e.g. 1391). Therefore we will first prove the Theorem in the manifold case in the C 2 
metric category, and then use perturbation estimates to obtain the result for g £ C a . To this 
end, we will often have dependence on \ \g\\ a even though we will be (for a specific Lemma or 
Proposition) assuming the g £ C 2 . 

Notation. 

• In what follows, we will write f{x) < C i,...,c„ g(x) if there exists a constant C depending only 
on Ci, . . . , c„, and not on /, g or x, such that f{x) < Cg(x) for all x (in a specified domain). 
We will write f(x) ~ Cl ,..., c „ g(x) if both f(x) < Cl ,..., c „ g(x) and g(x) < ci ,...,c„ f(x). If f,g 
take values in R d the inequalities are intended componentwise. We will write a b if 
Cib < a < Cib (componentwise for a, b vectors). 

• In what follows we will write d p K t {-, •) to denote the partial derivative with respect to the 
second variable of a heat kernel at time t. 



3.1 The Case of g e C 2 . 

We note that even though we assume g £ C 2 , we only use the C a norm of g. The idea of the proof 
of Theorems 12.1.11 and 12.2.11 is as follows. We start by fixing a direction p\ at z. We would like 
to find an eigenfunction ipi 1 such that ISp^iJ > Rj 1 on B Ci r z (z). In order to achieve this, we 
start by showing that the heat kernel has large gradient in an annulus of inner and outer radius 
~ -Rj 1 around y\ (j/i chosen such that z is in this annulus, in direction pi). We then show that the 
heat kernel and its gradient can be approximated on this annulus by the partial sum of (|2.1.ip over 

eigenfunctions ip\ which satisfy both A ~ R~ 2 and R z 2 \\ 1 P\\\l 2 (b ciRz {z)) ^ 1- By the pigeon-hole 
principle, at least one such eigenfunction, let it be tfi 1 has a large partial derivative in the direction 
P\. We then consider and pick pi _L and by induction we select fi 1 , . . . , ipi d , making sure 

that at each stage we can find tpi k , not previously chosen, satisfying \d Pk ipi k \ ~ Rj 1 on B Ci r z {z). 
We finally show that for the proper choice of constants 71, ..,7d < 1, the map $ := (71 y^, ..., r Yd ( Pi d ) 
satisfies the desired properties. 

When working on a manifold, we assume in what follows that we fix a local chart containing 
Br z (z), as at the beginning of section |2~21 

Step 1. Estimates on the heat kernel and its gradient. Let K be the Dirichlet or Neumann 
heat kernel on fl or M. , corresponding to one of the Laplacian operators considered above associated 
with g and the fixed a. 

Assumption A.l. Assume g E C 2 , and let a £ (0, 1] be given and fixed. Let constants 60,61 > 
depend on d, c m i n , c max , ||g|| Q , ex. We consider z, w £ Q satisfying 4f i? z < t? < 6iR z and \z — w\ < 
6 a R z . 

Remark 3.1.1. Provosition \ 3.1.~B below makes no assumptions on the finiteness of the volume of 
M and the existence of C coun t ■ It is also used in the proof of Theorem \2.2.8i 

Proposition 3.1.2. Assume Assumption A.l, 60 sufficiently small, and 61 is sufficiently small 
depending on 6q. Then there are constants Ci,C2,C' 1 ,C' 2 ,Cg > 0, that depend on d, 5q, 61, c min , 
Cmax; -Rzllfflla? ct, such that the following hold: 

(i) the heat kernel satisfies 



K t (z,w) ~g ; 



(3.1.1) 
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(ii) if ^SqR z < \z — w\, p is a unit vector in the direction of z — w, and q is arbitrary unit vector, 
then 

\VK t (z,w)\ ~g and \d p K t (z,w)\ ~gj t^ ^ (3.1.2) 



t 



d q K t {z, w) - d q Kf (z, w) 



<C 9 t^^- 
t 



(3.1.3) 



where Cg — > as Si — > (toit/i <5o fixed); here, Kf (z, w) is the usual Euclidean heat kernel. 
(Hi) if ^SqR z < \z — w\, and q is as above, then for s < t, 



K s {z,w) <c 2 t 2 
(iv) Ci, C2 both tend to a single function of d, c 



VK.(z,w)\ < c ^ t^?f and \d q K s (z,w)\ < c , t^?f; (3.1.4) 



min ; Cmax j Hslla, as tends to wii/i <5 fixed; 



The reason for the factor of R" which we have in i?"||<?|| Q above is to get scaling invariance. Propo- 
sition 13.1.21 is proved in subsection 13.3.11 for the Euclidean case and in subsection 13.5.31 

We continue with the proof of Theorem 12. 1.11 and !2.2.1t From here on, unless explicitly stated, 
we assume the existence C coun t. We have the spectral expansion 



+00 



K t {x,y) =J2 e A, Vj(aOv»j(l/)- 
3=0 



(3.1.5) 



Remark 3.1.3. The assumptions of Theorems 1 2. 1.1\ and \2.2~A\ say that \M\ — 1 (manifold case) 
or |fi| = 1 (Euclidean domain case). Thus, unless explicitly stated, we will assume in the lemmata 
below that we have R z ^d.c min .c max ,||g|| Q ,a 1- 

The following steps aim at replacing appropriately chosen heat kernels by a set of eigenfunctions, 
by extracting the "leading terms" in their spectral expansion. 

Step 2. Heat kernel and eigenfunctions. We start by restricting our attention to eigen- 
functions which do not have too high frequency. Let 



A L (A) = {Aj- : X 3 < At' 1 } and K H (A') = {X 3 : Xj > A't' 1 } = A L (A') C 



(3.1.6) 



A first connection between the heat kernel and eigenfunctions is given by the following truncation 
Lemma, which is proved in subsection 13.4.21 

Lemma 3.1.4. Under Assumption A.l, for A > 1 large enough and A' < 1 small enough, depending 
on Sq, Si, Cj., C'2, C' ly C' 2 (as in Proposition \3.1.2\) , there exist constants C-i,C^ (depending on A, A' 
as well as d, c m ; n , c max , ||<?||aj ol) such that: 



(i) The heat kernel is approximated by the truncated expansion 

K t (z,w)~c 4 3 ^j( z VjH e ~ 

j£A L (A) 



(ii) If ^SqR z < \z — w\, and p is a unit vector parallel to z — w, then 

d p K t (z,w) Y fj(z)d p ip j (w)e- 

jEA L (A)DA H (A') 



Furthermore, 



^(z)Vip J (w)e-^ t 



< C w t~ 



^Rz 



(3.1.7) 



(3.1.8) 



(3.1.9) 



where C10 — > as A — ► 00 and A' 



0. 
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(Hi) C3,Ci both tend to 1 as A — > oo and A' — > 0. 

This Lemma implies that in the heat kernel expansion we do not need to consider eigenfunctions 
corresponding to eigenvalues larger than At^ 1 . However, in our search for eigenfunctions with the 
desired properties, we need to restrict our attention further, by discarding eigenfunctions that have 
too small a gradient around z. Let 



A E (p, z,R z ,5 ,c ) 



\j e cr(A) 



1 

CO 



Rg\d p ipj(z)\ > 



B(z,±5 Rz) 



<p 3 {z'ydz' 



(3.1.10) 



Here and in what follows, f A f = \A\ 1 f. f. The truncation Lemma \3 . 1 . 41 can be strengthened, 
on average, into 

Lemma 3.1.5. Assume Assumption A.l, 5q sufficiently small, and Si is sufficiently small depending 
on 5q. For C%,Ci close enough to 1 (as in Lemma \3.1.4\ ), and cq small enough (depending on 
C2,C' 1 ,5o,Si) there exist constants C§,Cq (depending only on C3, C4, Cg, and cq) such that if 
■^5qR z < \z — w\, and p is a unit vector parallel to z — w, then 



\d p K t ( w ,z)\ ~g; 



ip j (w)d p (p j (z)e 



-A,t 



A J eA I ,(A)nA f f( J 4')nAE(z,^,(5o,co) 



(3.1.11) 



Step 3. Choosing appropriate eigenfunctions. 

Define the constants r y Vj as 



1 Vi = f w(zTdz'\ . (3.1.12) 

Note that since tfj is L? normalized, we have j Vj > R d J 2 . A subset of these constants and cor- 
responding eigenfunctions will soon be chosen to give us the constants {7^} and corresponding 
eigenfunctions {<PiA in the statement of Theorem 12.1.11 and Theorem 1 2. 2. 11 

The set of eigenfunctions with eigenvalues in Al{A) D Ah (A') n Ae{p, z, R z , So, c ) is well-suited 
for our purposes, in view of the following: 

Lemma 3.1.6. Under Assumption A.l, for 5 small enough, there exists a constant C7 depending 
on Co and C% depending on {So,c mm ,c max ,\\g\\ a ,a} and a constant b > which depends on Cq, d, 
Cmin; c max; ||s|| Q , a such that the following holds. Letp be a direction. For all j £ Ae(p, z, R z ,5q,Cq) 
we have that for all z' such that \ \z — z'\\ < bSoR z 




\d P vAz')\~c a Rz\i tf\ ■ (3.1.13) 



Moreover, there exists a index j e Al(A) D Ah (A') (1 Ae{p, z, R z , 5q, cq), so that we have 

l Vj < (C count )i (3.1.14) 
with constants depending on A, C±, C[, C2, C 2 , Ce, {d, Cmi n , Cmaxi ||ff||a:) &}> <5o, <5i . 
We can now complete the proof of Theorems 12.1.11 12.2.11 
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Proof of Theorems 1 2. 1 . 1\ and \2.2.1\ for the case g £ C 2 . Lemma 13.1.61 yields an eigenfunction that 
serves our purpose in a given direction. To complete the proof of the Theorems, we need to cover 
d linearly independent directions. Pick an arbitrary direction p\. By Lemma 13.1.61 we can find 
jx S A L (A)nA H (A')nA E (p, z, R z ,Sq,c ), (in particular ji ~ t^ 1 ) such that | j Vji d Pl ip^ (z) \ > cqR^ 1 . 
Let P2 be a direction orthogonal to (z). We apply again Lemma f3. 1 .6|. and find j2 < At^ 1 so 
that \"/ip j2 d P2 (pj 2 (z)\ > cqR^ 1 . Note that necessarily j2 ^ ji and ^2 is linearly independent of pi . In 
fact, by choice of p%, 

d P2 ip n =0. 

We proceed in this fashion. By induction, once we have chosen ji, . . . ,jk (k < d), and the corre- 



sponding pi, . . . ,pk, such that 



> cqR z , for / = 1, . . . , k, we pick Pk+i orthogonal 



> 



to ({*V<pj n (^)}n=i,...,fc) and apply Lemma T3.1. 6[ that yields jk+i such that J(p jk 1 d Ph+1 (pj k+x (z) 

CqRz 1 . 

^From here on we denote by 7* = 7^. for simplicity of notation. These are the constants {7;} 
appearing in the statement of Theorem 12.1.11 and Theorem 12.2.11 
We claim that the matrix 

A k +i := (lmdp n ( Pj m (z)) m ,n=l,...,k+l 

is lower triangular and {pi, . . . ,pk+i} is linearly independent. Lower-triangularity of the matrix 

follows by induction and the choice of Pk+i- Assume Yln=x a nPn = 0, then (X)n=i a nPm li^fji ( z )) — 
for alH = 1, . . . , k + 1, i.e. a e M fc+1 solves the linear system 

A k+1 a = 0. 

But Ak+i is lower triangular with all diagonal entries non-zero, hence a = 0. 
For I < k we have (V^-, (z),pfc+i) = and, by Lemma 13.1.61 

|(7;V^,pz)| >i?J x . 

Now let $fc = (71 (pj! , . . . , Jk<Pj k ) and $ = <&d- We start by showing that 

\m\ z (w-z)\\> d -^\\w~z\\. 



Indeed, suppose that 



|V$ fc |»(u;-*)|| < ^-\\w-z\ 



for all k = 1, . . . , d. For c small enough, this will lead to a contradiction. Let w ~ z = a/p;. We 
have (using say Lemma |3 . 1 .6[) 

\\^ k \ z (w -z)\\ = ||5>0ft**l*H = ||$>a K <I> fe L|| > (|a fc | - c£ |o,| J ^- . 

I l<k \ Kk / z 

By induction, \a k \ < X)f=i c 'll u ' — z \\- For c small enough, |cij| < U^-ili, This is a contradiction 
since || J2i a iPi\\ = ll w _ z ll an d | bi|| = 1< 
We also have, by Proposition 13.4.1] 
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Finally, by ensuring "'H is smaller then a universal constant for i = 1,2, we get from equation 



(13. 1.151) 



||$K)-$(u; 2 )|| = 



V$| fU)1+ (i_t) W2 (tyi - w 2 )dt 



ttUi + (l — t)t02 



Vi>U)) Oi -w 2 )dt 
f l 1 1 

~ / S-\\ W 1 ~ W 2\\dt > —CoWlV! - w 2 \\ , 

JO ^2 



which proves the lower bound (|2.1.5p . To prove the upper bound of (|2.1.5[) . we observe that from 
Proposition 13 . 4 . ll we have the upper bound 

l7^ Pl ^(z)|<i?- 1 

This completes the proof for the Euclidean case. 

We now turn to the manifold case. Let R z be as in the Theorem. We take c\ < Mo chosen so 
that 

\g a (x) - S a \ = \g u (x) - g a {z}\ < \\g\\ a \\x - z\\ a < e (3.1.16) 

for all x G B 2ci r z (z). For this g, the above is carried on in local coordinates. It is then left to prove 
that the Euclidean distance in the range of the coordinate map is equivalent to the geodesic distance 
on the manifold. We have for all x,y £ B Ci r z (z) 



dM(x,y) < 



f 


x-y 


dt < f 

9 JO 


x — y 


Jo 


\\ x - y\\ 


\\ x ~ v\\ 



(l + \\g\\ a t a )dt< a (l + \\g\\ a )\\x-y\\ 



The converse can be proved as follows. Let £ : [0,1] M be the geodesic from x to y. £ is 
contained in B 2c t M r Xty \(x) on the manifold, whose image in the local coordinate chart is contained 
m B 2 ( 1+ \\g\\ a )d M (x, v )(x). We have 



d M (x,y) = y IK(*)|| tf >(l-|Mla) jm)\\^>(l-\\9\\ a )\\x-y\\ 



□ 



3.2 The Case of g e C a . 



Proof of Theorems \2.1.1\ and \2.2A\ for the case g £ C a . We can now give a short proof for the g G C c 
case, relying on the C 2 case. We need the following Lemma, which we prove in section f3. 5. 21 

Lemma 3.2.1. Let J > be given. If 

\\g l n - g ll \\L^{B R {z)) ->7i o 

with ||c q uniformly bounded and with fixed ellipticity constants (as in (|2.2.1jl ). then for j < J 



\\<Pj ~ W,n\\L°°(.B R (z)) ->n , 
|| V((pj - <Pj,n)\\L°°(Bn(*)) °' 



(3.2.1) 
(3.2.2) 



|Aj — Xjn\ 



0. 



(3.2.3) 



Now, to conclude the proof of the Theorem for the C a case, let J = c${d, ^c m i n , 2c max , | \g\ \ a ,a) ■ 
R~ 2 . We may approximate g in C a norm arbitrarily well by a C 2 (A4) metric. By the above Lemma, 
and the Theorem for the case of C 2 metric, we obtain the Theorem for the C a case. □ 
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3.3 Heat kernels estimates 

This section makes no assumptions on the finiteness of the volume of M and the existence of C coun t ■ 
3.3.1 Euclidean Dirichlet heat kernel estimates 

We will start by proving the heat kernel estimates of Proposition 13.1.21 for the Dirichlet kernel K Q . 
These estimates are in fact well known, and we include their proof here for completeness, and also to 
introduce in a simple setting the kind of probabilistic approach that will be used to obtain estimates 
in a more general context. 

Proof of Proposition \3.1.2\ for the Euclidean Dirichlet heat kernel. Let below be a Brownian path 
started at point zell, and t(uj) its first hitting time of <9f2. We recall that as a consequence of the 
Markov property we have (see e.g. [22], eqn. (9.5) page 590) 

K?(z,w) = Kf(z,w) - E w (Kflr M (BZ(T(u)),w)xt>«u)) (3-3.1) 

Then, 

E w (x t Rd rM (i?^T(c i ;)(a;)), W )xt> rMM ) = E w ((4n)^(t - r{u))^e^-^ Xt>r( w )) 

/. .-d,.-d -||b(t(^))-^|| 2 \ / d ..-d -(i-ap) 2 « 2 

< E w I (4tt) 2 (i) 2 e n Xt>T(u,)j < ^ I (4tt) 2 (t) 2 e « Xt>r( w ) 



-(l-ao)-'fi 



-(i-s ) 2 fi 2 _ rf , -(i-gp) 2 



< (4n)~(t)~e « < (47r)~(i)~e ^i^) 2 = (4tt)— (£)~ 



(3.3.2) 

where for the first inequality we require ^ B ^~ w ^ sufficiently large, which is implied by choosing 
Sq < 1 and Si small enough. The last term can be made arbitrarily small by choosing 8\ small 
enough, independently of S (as long as, say, S < |). We also have 

Kf d (z,w) = (4?r)^ '(t)Tr e^T?^- < (47r)^(i)^ 

and 

_-2 

T ,B d . \ ,. . -i . . -d -|l*-™ll 2 ., .-d,.-,l - s o R 1 ,, . -d . . -d — # o Q\ 

= (4vr) — e — n — > (4tt) — (<)~ e — n— > (4tt) — (t) — e *? . ^ d 

If So < |, then ^ 1 ~^°- > > 5q and so by reducing S± (while fixing Sq) we can make the left-hand side 
of (|3.3.2|) arbitrarily small in comparison with (|3.3.3|) . Now, from equation (|3.3.ip we get (|3.1.ip for 
the Dirichlet kernel. 

-d -d -|I-BJ(t(^))II 2 

Note that the range we have for t and \\B*(t(u)) — w\\ imply that (47r)~ (t)~ e « is 
monotone increasing in t. Hence we also have 

(x 8 > T ( u )^ T(u) K(r(^)),»)) < {A^{t)^e ZS W- . (3-3.4) 
If we also have |& is large enough, then Kf (z,w) is monotone increasing in s, and therefore 

Kf\z,w) < Kf\z,w) < (4tt)tV^ . 

For and fixed So, we may reduce Si so that by (|3.3.4|) is small, and thus we obtain the first part of 
(f3~l~4f from [pIXTJl . 
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We now turn to estimates (|3.1.2|) and second and third parts of (|3.1.4[) . We differentiate equation 
(|3.3.ip and then we bound as follows: 



\VA\Xt>r( u )^ M mr(u),v3))]\\ 



, , U -J -||B*(T(q,))-„ 

Xt>r(c)(47r)— (t-T(w))-r e 



< 



< 



, X ^ , , „-<! -||B*(t(u,))-„ 

Xt>T( w ) v »( 4 ' r )~ r ( t - r(w))— e *-^» 



Xt>T(0.) oA =7^ (47T)-»-(t - T(w))T- e 4 (t - TM) 



Xt>r(w) " 

Xf>r(w) _ 

, — d -d-1 

= (4tt) — t _ 5-E t , 

, -d -d-1 

< (4tt) — 5-E a 



2(t-T( W )) 

|55(r)-HI 



(47r)~(t - r(w))~e 4 <*— > 



2(t-r(w)) 



2f 



(47r) 2 t 2 e 4* 



Xt>r(^) — e 



X*>t(w)" 



2 V* 



26i 



4«f 



25i 



c(<y ,«i) 



(3.3.5) 



where for the second equality we use the dominated convergence Theorem, for the inequalities in the 
fifth and in the penultimate line we choose So < 1 and S\ small enough. Note that 5\ — > implies 
that C(5o, Si) — ► 0. Observe that these estimates hold also with Vu, replaced by -j^. 
We also have 

V^iff (z,w) = V u ,(4tt) — t~ e — « — = (4tt) — -=^e — « — 

-d -d-l &)^?z / x 

> (4tt) — t— ^-i e T ( z - iu) 

(with inequality understood entrywise) where as above the last inequality holds for Sq < 1 and 5i 
small enough. If p is parallel to z — w, the same estimates hold if we replace \7 W by d p . Hence, for 
any fixed So, by reducing Si we get 

E> -(l-^o) 2 

(1/.. - \ t-7 7-^K d /_ ...\n ^ / a _\^-"* ^ _ IST 1 — 



and therefore 



|V^ 2 (z, W )- V m if* (z,w)\\ < (47r)T-_£ t T e 



The estimate (|3.1.2[) involving <9 P is proven analogously. The second and third parts of (|3.1.4p follow 
as above. Finally, to prove (|3.1.3p we use (|3.3.5[) to obtain 

\\d q K?(z,w)-d q Kf (z,w)\\ < (47r)^^t^e _ ^ _ < C 9 (S ,Si)t^^ . 

□ 
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3.3.2 Local and global heat kernels 

In this section, let K be the heat kernel, Dirichlet or Neumann, for: 

(i) a domain f2 (possibly unbounded), and a uniformly elliptic operator A as in (|2.2.3p . with 
geC 2 (Cl); 

(ii) a manifold Ai with g 6 C 2 satisfying the requirements in section l2"T2l and let A be the associated 
Laplacian. 

Remark 3.3.1. We emphasize that in this section we do not assume that the volume of M. is finite. 

Observe that in both settings the existence of an associated Brownian motion is guaranteed ( [39] 
for the M. d case and the manifold case then follows from uniqueness). The following result connects 
K with the Dirichlet kernel on a ball, associated with A, to which the estimates of the previous 
section apply: this will allow us to extend estimates for the Dirichlet heat kernel on a ball to the 
general heat kernel K . A more detailed account of the ideas in the following proposition appears in 
Stroock's recent book [50] (section 5.2 Duhamel's Formula). 

Proposition 3.3.2. Let z £ Q, and r < dist(z, dil), or z <E M. and r < rjj(z). Let x,y € B{z, ^r). 
For each path ( starting atx), we define t\ (lj) < ti (u>) < . . . as follows. Let t\ (lo) be the first time 
that B% re-enters B(x, |r) after having exited B(x, ^r) (if this does not happen, let ti(u>) = +oo). 
Let xi = B*(ti). By induction, for n > 1 let T n (uj) be the first time after t„_i(w) that B^ re-enters 
B(x,^r) after having exited B(x,^r), or +oo otherwise. Let x n — _BJ(r„). If T n {uj) = +oo, let 
T n+ k(oj) = +oo for all k > 0. Then 



K s (x,y) = K°(x,y) + Y,^ Cm^M.K) t ™ < s P ( r ™ < s ) ■ (3-3-6) 



n=l 

where 



Kf = K s 

is the heat kernel at time s for the ball B(x, hr) with Dirichlet boundary conditions. Moreover there 
exists an M = M(c m i n ,c max ) such that 

r 2 

P{r n < s) <d, Cmin , Cmax exp{-n— } . (3.3.7) 

M s 



Remark 3.3.3. In our applications of this proposition, we have r ~ S Q R Z . In that case, if si < 
5\R Z , for |& sufficiently large (depending only, on d, c m i n , c max ), the factor exp{— n-p^} can be made 
arbitrarily small. This gives us control (exponential in n) on the right-hand side of (|3.3.7[) . Hence, 
for any fixed Sq, for Si — > the right-hand- side of equation (|3.3.6[) is dominated by the first summand. 



Proof of Proposition VJ.S.'A The proofs for the case of a domain fl and the case of a manifold M. are 
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identical. We have, for any fixed, small enough, e > 0, 

P Q (BZ(s) E B e (y)) = P n (BZ(s) G B e (y), n > a) + P n (BZ(s) e B e (y),n < a) 
K°{x,y')dy' 

P u (BZ(s) e B c {y),r 2 > s\ n < s)P( n < a) 

P n (BZ 1 (s)eB e (y),r 1 <s,r 2 < S ) 
Kf{x,y')dy' 

+ P n {B x J{s ~T 1 )eB e (y),r 2 > S \r 1 < s) P( n < a) 
P n (B^(s)eB e (y),r 1 < S ,T 2 <s) 

K?(x,y')dy' 



Be(y) 



+ f E u [K°_ T1 (xi,y')\n < s] dyP(n < s) 
JbAv) 

+ P n (B^(s)eB e (y),T 2 < 3 ) 



Kf{x,y')dy' 

B4v) 

+ V / E w [K°_ T (xi,y')\T n < s] dy' P(r n < s). 

By dividing by \B(y, e)| and taking the limit as e — > + , we obtain (|3.3.6[) . 
In order to estimate P(r„ < s), we need the following 

Lemma 3.3.4. Let CI be a domain corresponding to a uniformly elliptic operator as in (]2.2.3[) . Let 

t be the first exit time from B R (z) C 51 for the corresponding stochastic process started at z. Then 
there exists M = M(d,c min ,c m&x ) > such that, 

P z ({r < a}) < d , Cmm , Cmax exp{-R 2 (2Ms)- 1 } . (3.3.8) 

Similarly for z G M and R < rjj(z). 

Proof. First note that without loss of generality we may replace by B 2R (z) with Dirichlct boundary 
conditions, and then replace B 2R {z) by M. d by extending the coefficients g^ to g 1 ^ defined on all of 
M. d . Let K be the associated heat kernel. For any s' > and x,y 6 B R (z) 

s'-i exp{ _^_^ } < Cmin!Cmaxid K sl (x,y) < d , Cmin , Cmax a'-i exp{-^-^} . (3.3.9) 
A\s' A 2 s 

holds for M = M( Cmin ; c max , d) and Ai — ^4j(c m i n , c max , d) (see |17| . Corollary 3.2.8 and Theorem 
3.3.4). We now follow a short proof by Stroock [ST]. By the strong Markov property, we have 

P n (B*Js) i B R {z)) = K [p(B^\s-r{u)) $ B r {z)) X{t{uj)<s} 

^From the lower bound in equation (|3.3.9p we have that if x £ dB R (z) and s > then P(B*(s) ^ 
Br(x)) > e(c m ; n , c max , d). Combining this with the upper bound in equation (|3.3.9[) we have 

P z {t <s)< £ (d,c wia ,c miai )- 1 P{B*{a) i B R {z)) < d ,c mln ,c max exp{-i? 2 (2Af s )- 1 } . 

□ 
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We go back to the proof of Proposition 13.3.21 To estimate P(t„ < s|r n _i < s, . . . ,n < s), we 
observe that between r„_i and r„, the path u> has to cross both dBz r (z) and dBi r (z): let t* and 
t** be the first time this happens, and let y* = Then 

P(T n <s\T n -i<s,...,n<s)<P(T**-T*<s)< sup P^'l sup | \B U («') - y*| | > ) 



J*eflj (z) \s'G[0,s] 



S mi „,C max ,d e 2A f3 

Therefore we have 



P(r„ < s) = p(n < s, r 2 < s, . . . ,r n < s) = ( Yl p (n < s \n-i < s,...,n < s) \ P(n < s) 



\l=2 / 

^c min , Cmax ,d exp{-n Qr) (2ms)- 1 } . 

(3.3.10) 

Renaming 128M to Af we get the lemma. □ 

Remark 3.3.5. As it is clear from the proof, the proposition holds for any boundary condition on 
a manifold or domain. 

3.3.3 Euclidean Neumann heat kernel estimates 

We use the results of the previous two sections to prove the Neumann case of Proposition 13.1.21 

Proof of Proposition \3.1.2\ for the Euclidean Neumann heat kernel . The starting point is Proposi- 
tion [33^21 which allows us to localize. We use Proposition 13.1.21 for the case of B^s rA z )- For this 
proof, we denote by C^-B] be the C2 constant for the Dirichlet ball case. For s < t we have using 
equation (|3.3.10[) . 



\K s (x,y) - K* s ° R *(x,y)\ = | [Kf_^ (x n (u), y)\r n < s] P w (t„ < s)\ 

n=l 

<c 2[B] f:t-UM-n( 6 ^f) 2 (Ms)^} 



(3.3.11) 



equation ]3.3. 101 
2 



<C 2[ B], So , Sl exp{- (^L\ (Ms)- 1 } 



This proves (|3.1.ip and the first part of (|3.1.4|) (see Remark l3.3.3[) . For the gradient estimates, i.e. 
(|3.1.2p . (|3.1.3[) . and the second and third part of (|3.1.4J) . 



T„ < S 



PjT n < S) 



\V y K s (x,y) - V v K?° R *(x,y)\\ < £ |v„l^ [k* s J*» (x n (u) , y)\ 

n=l 

**[*] £ t-i^L ex P {-n (^) 2 {Ms)--} 



n—1 s 

equation ]3.3. 101 

<C^B],6 ,5 x t * — exp{l-— ) ( MS ) > 
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giving us C9. By Remark 13.3.31 the exponential term from equation (|3.3.10l) can be made small 
enough so that we obtain estimate (|3.1.2|) as well as the second and third parts of (|3.1.4p . □ 

The proof for the manifold case is postponed to Section 13.5.31 



3.4 Heat kernel and eigenfunctions 
3.4.1 Bounds on Eigenfunctions 

We record some inequalities that will be used in what follows. 

Proposition 3.4.1. Assume g S C a . There existsbi < 1, andCp > that depends ond,c m i n ,c max ^ \\g\\ a ,ot 
such that for any eigenfunction ipj of Am on Br{z), corresponding to the eigenvalue Xj, and R < R z , 
the following estimates hold. For w € Bb 1 n(z) and x,y € Bb 1 n(z), 

l^Hl^CpP^R 2 )!^ {tpA (3.4.1) 



||V y -,(»|| : ■■ l . '' A )f : U \<pA (3.4.2) 



\\V^{x) - VipMW < Cp ^wP- (-f l^l 2 ) Ik-tfir (3-4.3) 



BrO) 



where P\{x) = (1 +x)i + P, P2(%) = (1 + x)% + P , P3 (a;) = (1 + x)^ + ^ , with (3 the smallest integer 
larger than or equal to ^j^. 

We postpone the proof to Section 13.5.21 Related estimates can be found in [T6J [48l H3 [46], [58] 
and references therein. 



3.4.2 Truncated heat kernel and selecting eigenfunctions 

The goal of this section is to prove Lemma 13.1.41 and 13.1.51 All the results of this section and their 
proofs will be independent on whether we talk about the Dirichlet or Neumann heat kernel, and on 
whether we talk about the standard Laplacian or about a uniformly elliptic operator satisfying our 
usual assumptions and whether we talk about a manifold M. or a domain Q. This is because the 
only tools we will need to obtain the results in this section are the spectral expansion of the heat 
kernel (|3.1.5j) . the elliptic estimates of Proposition 13.4.1] the assumption on C 'count (|2.1.2|) . and the 
bound 

K t (z,w) < Kt(z,z)*K t (w,v))$ . (3.4.4) 
which is a straightforward application of Cauchy-Schwartz inequality to (|3.1.5p . 

Proof of Lemma \3-l-4\ We upper bound the tail of the heat kernel: 

(3.4.5) 

by (|3.4.4j) and Proposition 13. 1.21 For A large enough, this implies (]3.1.7|) . 



<e~ \Kdz,z)\*\Kdw,w)Y- 
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For the gradient, we use Proposition 13.4. l\ and observe that x n e * is a decreasing function if x 
is large enough. We let r w = (1 — Sq)R z . 



\ i >At- 1 



< 



e 3 E V>j(^V>j(w)e' 



i «;)ilV 



it- 1 



e 2 



2(Ai*) 2 -e 



1 - A .7 



< Cp , 5ie ^|^(z,z)|M E / 



. -A -d 1 



Kt(w' , w')dw' - 



Now we consider the contribution of the low frequency eigenfunctions to the gradient. Proceeding 



as above, and recalling that in this case A,r^ < 4-, we obtain 



-A,t 



2 ^2 -2 

r r 



v Aj<A't-! 



JB(w,±r w ) 



2 e -A,t 



< Cp |^ t (z,z)|5- [ 2 / |^| 2 P 2 (A,-^) 2 e-^* 

7-10 \A J <A't-l </B Kj r ") 

<opfl.f^)|ir*(*,*)|i-[ E / 



\ifij\ 2 e- 



■ [^) A.h::.:;^-| E f I-!"' 



<c P ,c 2 P2 [ ^f) —\K t (z,z)\i I / if t/2 X'(t»',to') 
, A! \ 1 d , i d 

<C P ,C 2 P2 U- — * • 

oi / r w 
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Thus, by reducing A' we get the bound (|3 . 1 . 8[) and (|3.1.9[) . (We note that an alternative approach 

Xj t 

to the introduction of e^^A 7 " would have been to reduce S± and note that ^-f- is as large as we want 



in comparison with — , and then to reduce A' to compensate for the reduction in 8\. 



□ 



For a domain with Dirichlet boundary conditions, we automatically have a bound on C COU nt as 
in (|2~L2l) : 

Lemma 3.4.2 (Weyl's Lemma for Dirichlet boundary conditions). Let Q. be a domain in M. d , A a 
uniformly elliptic operator on Q as in (|2.2.1|) . with Dirichlet boundary conditions. Let Ao < Ax < • • • 
be the eigenvalues of A. Then 



#{j : \j < A} <d, Cmm , Cmax MX* 



(3.4.6) 



Proof. Let K n the associated heat kernel. Extend the coefficient g % i to R d \ f2 by letting gr y = S 1 ^ , 
and let K be the associated heat kernel on M. d . Then K is pointwise dominated by K. Then by 
estimate (|3.3.9[) we have, following [23]: 



#{j : Aj < A} < e • > e~ = e • / if" (x, sc)da: <e- Ki (x,x)dx <d, Cmin , Cmax |0|A* 

7n A Jn x 



□ 



Remark 3.4.3. Notice that in the Dirichlet case C CO unt is independent of Q. In the Neumann 
case, if one has good enough estimates on the trace of the corresponding heat kernel, the same proof 
applies. In general these estimates will depend on £1, and C 'count will not be independent of fl (see 
e.g. mW$)- 



Proof of Lemma \3.1.5[ In view of Lemma 13.1.41 we will show that the terms in the eigenfunction 
series corresponding to j S AL(A)nAjj(.A') but j Ae(p, z, R z , So, cq) do not contribute significantly 
to the lefthand side of (|3.1.1ip . Let Ai = A L (A) n A H (A') n (A E (p, z, R z , S Q , c )) c . We thus have 



ip J (w)d p ip J (z)e A ^ 
l^iHI \d P ^(z)e- x ^ 



< 



< 




w)\ 2 e^ x ' t 



\°pw( 

,A 3 €Ai 



2l| 2 e~ A,t 



z)\ 2 e~ Xit 



< K t {w,w)* ( J2 
Aj-eAi 

/ 

\ A ' eAl S(z,i(5 it*) 

< c ^ K t (w,w)? -j K t (z',z')$ dz',. 



< K t (w,w)2 Co- 
il 



\ Vj \*e-™ 



Hence by reducing Co and using Proposition 13.1. Al together with Lemma l3.1.4l we conclude the proof. 
We note that the constant Cg comes into play since we have to estimate the left hand side of equation 
(13.1.111) . □ 
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Remark 3.4.4. The following proof is the only place where we use the bound on C coun t (|2.1.2[) . 

Proof of Lemma \3.1.6\ For sufficiently small b, Equations (|3.4.2[) and (|3.4.3[) from Proposition ^. 4.11 
together with the definition of Ag(p, z, R z , Sq, Cq) give equation (|3.1.13[) . We turn to showing equation 
(|3.1.14p . which is where C count will appear. 

Since at this point all the constants are fixed, to ease the notation we let A := Al(A) D Ah {A') n 
Ae(p, z, R Zl So, Co). Let w € B(z, 5qR z ) n B(z, ^8qR z ) with w — z in the direciton of p. Observe that 
Proposition 13.1.21 and Lemma [3.1.51 imply that 

K*(w,w)^- ~% ~g \d p K t (w,z)\ 



E iPj(w)d p tp J (z)e 
A 3 eA 



-Xit 



$Hc min ,c max ,|| S || ct ,a} Rz 1 Y] \^j( w ) e ~ 

Aj eA ' \Jb(z±s r») ) (3.4.7) 



A '' ' Wj? 



<R^[Y^H 2 e- 2X A E/ 

Va^ea / Va.ga^M^) 



<R- x K M {w,w)* [ / 

Va^a^M^) 



<Pj\ 2 



giving, with constant depending on Ci,C' 1 ,C 2 ,C' 2 ,C 6 ,{d,c min ,c max ,\\g\\ a ,a}, 



^(f)-^{ (3-4.8) 



Thus, by the pigeon-hole principle and Weyl's bound (|2.1.2j) . we have Aj S Ai with 



i /i?r 2 



Ccount \ t J JB ls _ (z) 



\<Pi\ ■ (3-4.9) 



This gives equation (|3.1.14p . □ 
3.5 Supplemental Lemmata for the Manifold case 

We will initially be interested in localizing the manifold Laplacian /Am to a ball B = Br(z), 
R < r u( z )i hi a coordinate chart about z, satisfying the assumptions in the Theorem. We will 
rescale up so that R — 1 (and rescale the volume of M. accordingly) . We impose Dirichlet boundary 
conditions on dB, and denote by A B this Laplacian, which has the expression (|2.2.3[) . We will 
compare A B with the Euclidean Laplacian A B on B (also with Dirichlet boundary conditions). 
We will then compare A B with the global Laplacian Am on the whole manifold (with Dirichlet 
or Neumann boundary conditions). The first comparison is most conveniently done through the 
associated Green functions. We use the following notation: 

(i) A B £j — fij^j is the eigen-decomposition of A B , with sorted eigenvalues < [Xq < fii < . . . , 
A B £j = is the analogous decomposition of A B , and A.M<Pj — Ajipj the one for Am- The 
eigenfunctions are assumed to be normalized in the corresponding natural L 2 spaces. 
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(ii) G B is the Green function on B, associated with A B , with Dirichlet boundary conditions, and 
K B the corresponding heat kernel; 

(iii) G B is the Green function on B, associated with A B , with Dirichlet boundary conditions, and 
K B the corresponding heat kernel; 

(iv) the quadratic form associated with g, restricted to B, will be abbreviated as 



g B (u,v) = 




(3.5.1) 



for suitable u, v. 

We will use estimates from [25], where they are stated only for the case of d > 3. Our Theorems 
are true also for the case d = 2 (and trivially , d = 1). This can be seen indirectly by considering 
A4 := M. x T and noting that the eigenfunctions of M. and the heat kernel of M. both factor. 

We let 

(L*) P {B) = {/ : B -> R measurable : < +00} 

where 

\\f\\ (L , )P =su V t\{xeB: \f(x)\ >t}\7. 

t>Q 

We recall the following Theorem from [35] 

Theorem 3.5.1. Suppose d > 3, and g £ L°° and uniformly elliptic with c mln and c max as in 
(|2.2.1|) . There exists a unique nonnegative function G B : B — > R U {00}, called the Green function, 
such that for each y £ B and any r > swc/i that B r (y) C 

G B (-,y)eiy c 1 ' 2 (i?\ J B r ( 2/ ))niy c 1 - 1 (i?), 

G B | ai3 = 0, and /or a// cj) £ C™(B) 

g B {G B {;y)^) = ^y). 

Moreover, for each y £ B , 

(i) G B (;y)£(L*)^, with \\G B \\ , < diCjnin l 

(L«)3=2 

(ii) VG s (-,y) G (£*)A with ||VG B || , <d,c mMa 1 

fmj G B (x,2/) > d , Cm „ iCmin \x - y\ 2 - d for \x - y\ < \d{y,dB) 
(iv) 

G B (x, y) < dlPm „, Cmin \x - y\ 2 - d (3.5.2) 
If g £ C a we also have (see page 333 in (25J) 
(v) 

V y G B (x, y) < d ,c ma „c mta> a,|| s || Q \x - y] 1 -* (3.5.3) 

(vi) 

\V x G B (xi,y)- V x G B (x 2 ,y)\ <d,c m ^,o aSn ,R«, a ,\\ g \U ^ _ ^p+^i + ^ _ y | d+a _i ■ ( 3 - 5 - 4 ) 
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Simple consequences of the bounds above are the following inequalities, which we record for 
future use: 

/ \G R (z,y)\ p dy < Cl ^d,c mm ,c^ P R {2 - d)+ * (3.5.5) 

Jc 1 R<\\z-y\\<c 2 R, 

\V y G R (z,y)\Pdy <cx ) c 2) d,c min ,c max) a,|| s |U,p , (3.5.6) 

ci-R<l|z-all<c2« 

\y y G R (z,y)\dy < d ^ Cin ^ Mla R, (3.5.7) 



which are an immediate consequence of f|3 . 5 -3[) , and are valid for c\ , C2 > and < R' z < R z . 

We recall that if we only assume uniform ellipticity, without any assumption on the modulus of 
continuity of g, then we have no pointwise estimates on VG. 

3.5.1 Perturbation of eigenfunctions 

We start by comparing eigenfunctions of the Euclidean A B with eigenfunctions of A B . We remind 
the reader that we have rescaled up to R = 1. 

Lemma 3.5.2. Let J > and n > be given. There is an £q = eo(J) so that if e < eo and 
Id : (B, <5 U ) — > (B, g l i) is 1 + e bi-Lipschitz, then for j < J, 

Uj - ij\\^(B) < V , - Mil < Wj 

Proof. This follows from Lemma 20 in [24]. □ 
Lemma 3.5.3. There is an integer [3i oc > such that the following bounds hold: 

and if g lJ G C a we also have 

||Vj,fj||z,°°(fl) S, Cmln ,c max ,a : || 5 || Q Hifiif 1 "" , {3.5.9) 

with [3i oc = for d odd and (3i oc = i for d even. 
Proof. By the definition of G, and by recalling that G\dB = 0, 

& = G B A B £j = fljG B £j = ... = 

k 

Let pi, qi be such that 

l<i<k 

Then, using Young's inequality we have 

|| G B * ■ ■ ■ * G B J \ L i k < qh _ imm \\ G B * ■ ■ ■ * G B J \ L i k -i\\G B \\ L P k 

k fc-1 

~q k -2,q k -l,Pk-l ■■■ i$?l,?2,P2 I|G B ||lpi •••||G'' B ||LPfc . 

We have G B G L 3 ^ by Theorem l3.5.1[ we take pi — 4^3 and take k — d — 1 and get 

<f d l x = (d- -(d-1) + 1 = 0, 
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for odd d 



and for even d 



Now, for odd d, 



d-ld-2 d-1 



1d- 



2 d-1 2 
dd-2 d d,d-2 



1 = - 
2 



2d- 1 2 



ll6l|oo<^- 2 | LG B *. . .*G B J |^|||f||^ <A,- a ||G B ||^ <A/ 



which gives the first desired bound. If d is even do the same with | replacing 
For the gradient estimate, we have 



d-1 
2 ' 



where we used the defining property of G B in Theorem 13.5.11 and Green's Theorem. We estimate 
the last term by (|3.5.8p and equation (|3.5.3p to get the desired result. □ 

We can now convert the L 2 -estimates in Lemma 13.5.21 into L°°-estimates. We will need the 
following 

Lemma 3.5.4. Assume that \g % ^(x) — <5 y | < e for x € B. Then for tp g C^°(B) we have 



(V (g b (z,w)-G b (z,w)) ,Vip(z))d 

and if g lJ € C a we also have 

jjy (g b {z,w)-G b {z,w)) Mi(z))dz 



Z/»(B) 



< ellV^II^UVVIU-^ (3.5.10) 



<_ 



as well as 



( V (G B (z, w) - G B (z, w) , V& (z))dz 



< 



L°°(B) 



rf,C m i n ,C n 



,«,|| fl |U e£/(w) ftoc (3-5.11) 



,a,||fl|U ew(w) ftoe (3-5.12) 



with fJioc as in Lemma \3.5.3\ . 
Proof. 



L°°(B) 



L°°(B) 



J (V (g b (z,w) - G s (z,w)) ,Vip{z))dz 

f V S ij di (g b (z, w) - G B (z, wj) d^{z)d 
Jb id 

(s ij diG B (z,w) - g' 3 (z)d 2 G B (z,w)^ d^{z)dz 
(g l3 {z)d l G B {z,w) - S ij diG B (z,w^j d J tP(z)dz 



L°°(B) 
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/ E - ^)d t G B {z,w)\ d^(z)dz 



h3 

< c||V tf G B ||£i||W|U-(B) 



L~(B) 



which gives (|3.5.10[) . Using Lemma [233] one also gets (|3.5.11[) and (|3. 5. 1 2[) . 



□ 



Lemma 3.5.5. Let J, r\ > be given. Let (3i oc be as in Lemma \3.5.3[ There is an eo which depends 
on J, 7], d, c max; c min , \\g\\ a , a), so that if e < e , and \g tl (x) -S ll \ < e, for x £ B, then for j < J, 



It 1 ] ~ Mil < Wj 



(3.5.13) 



110 - OlU°°(B) ^d, Cm i n ,c max ,||g|| ct ,c VQlifa) 

where Q± is a polynomial of degree 2/3; oc . If g £ C a we also have 

_ li)IU°° ~d,c min ,c alilx ,\\g\\ a ,a VQ^iP-l) 

where Q2 is a polynomial of degree 2(3i oc + 1 . 



(3.5.14) 



(3.5.15) 



Proof. The bound (|3.5.13[) follows from Lemma 13.5.21 Let and pt be as in the proof of Lemma 
We have using the definitions of G B and G B 

Zt(w)-£t(w) = / ^S i *a i G B (z,w)d j £ l {z)-g i *(z)aiG B (z,w)d j i l {z)dz 



1,3 



f (d i G B (z,w)d j ^)-d i G B (z,w)d^ l ( 

JB i,3 

+ (s ij diG B (z, w)d j l{z) - g^(z)d l G B (z, w)d^i(z)) dz 
£» + [ X/ (®iG B (z,w)djt;i( 

J B ■ ■ 



z)-g^(z)d l G B (z,w)d^ l (z)\dz 



= E 1 (w)+ (M -G B *ti(w)-fn-G B *tii(w) 

= E\w) + {m - fit) ■ G B * + fit ■ G B * (6 - lt)(w) 

= E^w) + E 2 {w) + fit ■ G B * (& - &)(w) . 



where we have from equation (|3.5.12|1 



^ 1 |U~(B) ~d,c mi „,c ma x,a, \\g\\ a ^f-lif-l) '" 



and 



\\E 2 \\ L °°(B) ^,d,c mia ,c m ^,a,\\g\U W ' (w) A ° c • 

Iterating, we have 

I6M-6HI = \E 1 {w) + E 2 (w) + fifG B *^i-^){w)\. 

= \E x (w) + E 2 (w) +fn-G B * (E 1 +E 2 +fn-G B * (f, - it)) (w)\ . 



(3.5.16) 
(3.5.17) 



A„c-i 



fc=0 



< \\E^+E 2 \\ L ^ {B) {vi\\G B \\mB)) +tf l "\ G B *...*G B X Si(w)-iiW)\ 

(Si oc times 
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k=0 



! \\l^b)) +tf loa \\G B \\ L ^Mi-iih 



< 2» ?W ( W ) /S,M E fil+n-it e '=vQi(fii) 



fc=0 



where we require for the penultimate inequality e < r/. 
To prove the gradient estimate, 

V(£ - &))(y)| = V, | £0«G*(z, j/) fl X&C0 -^^^(^jf'^^z) 

X7 y G R (z,y)^ 3 (z) - V y G B (z,y)^ j (z) 

V y {G R -G B )(z,y)\ ■ \fiMz) 
+ \V y G B (z,y)\ ■ fljijiz) - Hj£j(z) 
Now using equation (|3.5.14[) , Lemma [3. 5. 21 and Theorem 13. 5. II we get equation (|3.5.15[) 
3.5.2 Bounds on Eigenfunctions 



< 



(3.5.18) 

(3.5.19) 

(3.5.20) 

(3.5.21) 
(3.5.22) 
□ 



The main goal of this section is to prove Proposition 13.4.11 We note that the inequalities (|3.4.ip . 
(|3.4.2[) . and (|3.4.3[) arc invariant under scalings of the metric, and so, once again, we assume in the 
proof of this Proposition and in all the Lemmata that R — 1. In this section all constants subsumed 
in < and > will in general depend on d, c m ; n , c max , | \g\ \ a , a. We will need the following result. 

Lemma 3.5.6 (Lemma 3.1 from [25 ). Suppose h is a bounded solution of A B h = in B. Then 
|Vft(a;)| <d,c mta ,c mM£> || 3 |Ua (l-distfojOr^lfclli-^). (3.5.23) 

Lemma 3.5.7. Assume that g G C a and A B h = on B = B\{z). Then for any r < 1 

||/l|U°°(B r (*)) ^d,Cnto,Om« C r\\h\\ L l(B 2r ( Z )) 

and 

||Vft||ioo( S (*)) ^<i,O mi n,C mlut ,||ff||a,Q I I k I U 2 (B 2 r (*) ) ■ 



Proof. Let r as above be given. Fix < <ii < <Z2 < 1. By the coarea formula [23), we have 



\h{x)\dH d - 1 {x)dt= / |WG r |da; 

Oi J {x£B r (z):G r (x)=t} J{x£B r (z):ai<G r (x)<a 2 } 



< \\h\\LZ(B Az ))\\\VG r \\\ LHB 
^ \\h\\mB r (z)) 
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by estimate (|3.5.6p . Hence there exists t* G [01,02] such that 



\h(x)\dH d -\x) < \\h\\ L2iBAz)) (a 2 - a^ 1 



Now, by A B harmonicity, 



i*wi = 



'{yeB r (z):G-(y)=f} 
' IIVG 



h(y)^-(z,y)dn d ~ 1 (y) 



t-({cr(»)=f}) / \h{y)\dH (y) 



~ INU 2 (,B r (z)) i$ai,o 2 ,<<,c I>ll „,c 



£ 2 (-Br(*)) 



Essentially the same proof holds if we replace z by w G i? r (z)) in the above estimates, giving the 
desired bound on h. 

In order to estimate the gradient, we use Lemma 13.5.61 which gives us 



□ 



HfflU.a CrlNk~(B,.(z)) • 

which implies the desired estimate. 
Lemma 3.5.8. Assume that g G C a . Let £j and ip k be as above. Then we have the estimate 

M 3 <Pk\\ L -j^ {B) <d,c m , n ,c max ,|| S || Q: a (((Ai+Afc))3 +p, j )(p. j ) f}l "\\<Pk\\L*{B)- (3.5.24) 



Proof. By the Sobolev embedding Theorem it is enough to prove that 

l|V(|^fc)IU 2 (B) ^ (((Ai+Afe))* + Mj)(£i) A °1l<^IU 2 (B) 
To this end, first note that we may write 



A B (^ k ) = Vk A B ^+^A B ^ k + £ g^dmdjli 



1,3=1 



(3.5.25) 



(3.5.26) 



and so A B (£-,■</?&) is defined as a function, and not just a distribution. Observe that ipi does not 
satisfy any particular boundary condition on dB, however since £j = on dB, integration by parts 
gives 



(ijipk , A B (£• ^ ) ) B = 9B ( V (I -ip k ) , V <p k )). 
Now, since g is a positive quadratic form, 

(V(fjVfc),V(|j93fc)> ffs = (OVpfc + VkV£j,tjV<Pk + <PkV$j)g B 



(3.5.27) 
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and therefore 

9B{(ljV¥k,ijVipk) < 9B(V(ijtpk), V(^fc)) + 2\9R(<PkVii,iiV<Pk)\ 
= (i j <Pk,A B (i m )) B + 2\g R ( l p k V^^ j Vi Pk )\ 

Z (Ai+AfcJII^IIoollVfcll^fBJ + IIV^HoollVfclU'CBjflB^VVfc.^VvJfc)' 

< (/ij + Afc)(/i i ) 2 ^"' t =||^ fc ||| 2(s) + /i J (/ij)^ oc ||^fc||i2 (s) 5' J B(| J V^ fc ,| 7 V^A ; )* 



giving 



Finally, 



gBiZjVipk^jVVk)* = H£jV<£ft||z,2( BrW ) < ((fij + X k )2 + fi )(jl j f l '"=\\ip k \\ L 2 (B) 



V(^v»*)<IVfc||V&| + |C-VVfc| 
gives equation (|3.5.25p . □ 

Proof of Proposition \3.4-l] We recall that we rescaled so that i2 = 1. Let ip — aj^j be a (finite) 
sum of (Euclidean) Dirichlct cigcnfunctions of B such that 

- < ip(x) < 2 and x € B R/2 (z) C B 

and 5]i h i < C 1 , Mi < C, 1 < j < N. One may obtain such a sequence by taking -0' € C°°(B) with 
< ip 1 < 1, ^'|s R/2 ( z ) = 1 and ip'\dB(z) — and then take ip to be a truncation of the eigenfunction 
expansion of ip' . Let %p — dj£j be the sum of the corresponding Dirichlet eigenfunctions for B 
with respect to A B . By lemma [5.5.51 and \g lk (x) — 6 lk \ < e (with e sufficiently small), we have, for 
x € B R/2 (z), 

\ < fa) < 3 



By Lemma [3.5.81 



05,11 2d < 2d < > lajl H&^ill 2d < 



'L^(B)) ~ " JII L^(B)) ~~ ^— ' 1 111 rJ "L^(B 
X)N((£i + A.,')* +M 1 )(A 1 ) Aoc ||^j||l2 (b) <c (Aj+l)'||(Pj||i2(B) 



(3.5.28) 



We are now ready to prove inequality (|3.4.1j) . Let ro = i?=l>ri>r2>--->-| — § ■ Write 
Pjlflroto on as ^'IsroW = « + « where 

« = G S (AV)=A,-G S ( W ) (3.5.29) 

since G s is the Green function for the Dirichlet problem on B ra (z). Hence A B u — 0. We use (see 
below) Lemma [3.5.71 in conjunction with the above decomposition, to show that ipj € L°°(B rao (z)). 
We will then (see below) get (|3.4.2[) from differentiating (|3.5.29|) and using Lemma [3.5.71 Initially, 
by (|3.5.28j) . Theorem [3311 (|3. 5.29ft and Young's inequality, with p = £^ and 1 < pi = d _g+ m 
(with < rji < 4 of our choice, implied by the estimates on the Green function in Theorem 13.5. ip . 
we have 

IMUpiCB^O)) ^ ^j\\<Pj\\LPo{B ro {z)) 
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giving, by Lemma \'S. 5. 71 (since p\ > po > 2), 

\\u\\ L ~>(B ri (z)) < \\u\\L2(B ra (z)) < (l + Aj)||^3||x,2( flroW) . 

Thus, we have 

\\<Pj\\LPi(B ri ( Z )) ^ (1 + Aj)lkjlU"o(B ro ( z )) < (Aj + l) 5 ||<^j|U 2 (B) • 

Let 1 < f»i = d _ 2 _ 4 ,-^|-. - — — (with < 77,; < 4 of our choice) and u$ = G ri (A Ti ipj). Similarly, we 
have 

\\vi\\^(B ri _ 1 (z)) < ^]\W]\\h' , i-i(B r ._ 1 (z)) 

and for Ui — tpj — Vi 

\\ui\\L°°(B ri (z)) ^ ||w||i2( Bri _ iW ) < (1 + Aj)||^|| L 2( Br . 4 _ iW) . 
Thus, we have by induction 

\\tPj\\Lm(B r .(z)) < (l + AjOII^IU^-iC-Bn-iW) ~ + 1 ) t+3 llwlU m (B) • 

Let /3 be the smallest integer larger or equal than 2=2, We may choose {r/i} so that = oo. 
This gives equation (|3.4.ip . 

In order to upper bound ||Vy>j|| we note that (recalling that ~ R ~ 1) 

||V^|U-(B r „) = ||VG r "(A^^)IU»(B^) < \ 3 \\G^\\ L1{Br J<p\\ L ~ {Brp) . 



We also note that we have 



IIV^IUoc^^) < \\(p\\ L °°(B rfi ) 



from Lemma T3. 5 .71 Thus combining the last two estimates, we have (|3.4.2p . 

Finally, we prove (|3.4.3[) . Let X € C°°(R) be a function so that < X < 1 ; x( s )|s<A"i = and 
x( s )|s>k 2 = 1- We define 77, a cutoff function, such that r/\ B ^ z i R ^ = 1 and t]\\ x \>±r = as follows. 
Define rj{x) = x(G(z, x)), and choose K\, K 2 above so that 77 has the desired cutoff radius. We get 
that 

A B ( V )(x) = A( X (G(z, x))) = d Xi 9 ij d XjX (G(z, x)) = £ d Xi (x'(G(z, x)))g^d Xj G(z,x) 
= X "(G(z,x)) ( ^2 d Xi G(z, x)g ij d X] G(z, x)) + x'(G(z, x) J A x G(z,x), 



ho 



where the second term in the last line is as A x G(z,x) is a distribution which equals on the 
support of x'{G(z,x)). By choice of X and Theorem 13.5. U this gives A(rj) < 1. 
Now, let x,y G B = B(z, jR). Let G = G B be the Green's function for B. 



\\V<p(x) - V^(y)|| = \\V(<pr,)(x) - V^XiOII 



< 



y (W 1 G(x,w)-W 1 G(y,w))A B (^r 1 )(w)dw 
\A B { vv ){w)\ L ~ {B) [ \\V 1 G(x,w)-V 1 G(y,w)\\dw. 



(3.5.30) 
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We have (using uniform ellipticity as well as Proposition 13.4.1]) 

|A S M)H| < \r}& B <p(w)\ + \<p(w)A B v \ + \Vfi\\V<p\ 

< (HA^IU + A|MU) \\<p\\ L °°(B) + ||V??||oo||V^Uoo (B) 

< (1 + A) ||^||z,«.(s) + ||V(p|| L = (B) 

< (1 + A) Pi(A)|M| i2(B) + AP 2 (A)||^|| i 2 (B) 

< ((1 + A)Pi(A) + AP 2 (A)) y\\ LKB ) 



(3.5.31) 



by using (p?X2"j) . Combining (|3.5.30j) with (|3.5.3ip and (|B3^|) we get 

\\V<p(x)-V<p(y)\\ < P 3 (X)M\ L 2 {B) \x-y\ a . 



(3.5.32) 
□ 



Proof of Lemma \3.2.1\ This follows from Lemma 20 in [24] together with Proposition ^. 4.11 we have 
Qi+a functions which are close in L 2 (B) . Hence, they are also close in L°°(B), i.e. equation (|3.2.1|) 
holds and so does (|3.2.2| . □ 

3.5.3 Heat kernel estimates 

This subsection makes no assumptions on the finiteness of the volume of M. and the existence of 
C count for the manifold M. It will however use these properties for a manifold ball. 

We fix a ball B — Br(x) for which we estimate the heat kernel K B by comparing it to K B . 
Suppose that is an orthonormal basis for L 2 (B) (with manifold measure). In this section all 
constants subsumed in <, > and ~ will in general depend on d, c m - ln , c maX) ||<7||a) 

Lemma 3.5.9. Let Ax > 1 and a sufficiently small t/q = r)o(Ai) > be given. Assume 6q is 
sufficiently small (depending on r/o, A%, as well as the usual d, c m i n , c max , a), and \g lk (x) — 

5 lk \ < €q. For y S Br{x) C ft, with \x — y\ 2 < t ~ R 2 < 1 in a similar fashion to Assumption A.l, 
we have 



^ A i 

// in addition we also have \x — y\ 2 ~ t then 
The constants in (|3.5.33[) go to 1 as t]q — > 0. 



(3.5.33) 



< 



Ai,d,c m ^,c mixl ,\\g\\ a ,a Vo ' —t 2 ■ (3.5.34) 



Proof. We apply Lemma [3X5] with J = #{j : /x* < Ai/t} < (4 1 ) 2 R d ~ A{ and with r) < r/ . Let 
e be as guaranteed by Lemma 13.5.51 Since £j's and £j's are L 2 -normalized, Lemma 13.5.31 and 13.5.51 
implies for fj-i < ^ 



< Ux) - ii(x)\Uy)\e-^ + My) ~ ii{y)\Mx)\e-^ + Ux)\Uy)\\e-^ 

< Qi(A 1 t~ 1 R 2 )n (My)\e-^ + Mx)\e-^*) + Ux^Uv)^^^ 
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Using Weyl's Lemma (Lemma 13.4. 2[) for the ball with Dirichlet boundary conditions (see Lemma 
15X2)1 . we have 



E &(*)&(y)e^ 4i - E 
/^<-r Mi<-r 



< A? A "(l + Ai)tjJ 



by the (Euclidean) estimates in the proof of Lemma \3 . 1 . 41 and since i? < 1. We obtain the desired 
estimate (|3.5.33p by taking 770 sufficiently small. Similarly, 

|6(z)V&(2/)e-^ - 6(a;)V|i(»)e-^| 

< - ^WllV^d/Jle-*' + |V&(tf) - V^(2/)||^(x)|e-^* + |6(z)||V&(i,)||e-'* t - e~^\ 
<r l ((A 1 t- 1 R 2 )^\^M\e-^ t + (Axt- 1 Ry^ +1 ^^ 

Thus, equation (|3.5.34|l also clearly follows by ry sufficiently small. □ 

Lemma 3.5.10. Let t]q > be given and assumed to be sufficiently small. Assume eo is sufficiently 
small (depending on r]o, as well as the usual d, c m i n , c max , ||<?||a, a), and \g lk (x) — 5 lk \ < eo- For 
V € Br (x) C ft with \x — y\ 2 < t ~ R 2 < 1 (in a similar fashion to Assumption A.l) and s < t, 



Kt{x,y) ~ 7)0 ,d,c mla ,c mai ,||s|| ot ,a K?(x,y) , 
Kf{x,y) <^,d,c m ,n,«w,|lffllo..« K s( x 'V) > 



and 



\\VK?(x,y)\\ < no , d ,c min ,c mMa ,a f( sir2 r 2 ' 
// m addition we have \x — y\ 2 ~ t t/ien 

Vk?(x,y)-VKf(x,y) 

The constants in (|3.5.35p go to 1 as 770 ~~ >• 0. 
Proof. We estimate the tail : 



(3.5.35) 
(3.5.36) 

(3.5.37) 
(3.5.38) 



— \Xit 



<e~* Al E €i(as)fi(y)e-* wt 

using ma 



This, combined with ()3.5.33[) . for A x large enough, gives (|3.5.35|) . From [TB] we also get (|3. 5. 361) . 
We also have 



-fliS 



[ V x G B {x lW ) E A B &M&(y)e- 
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f v x g b (x : w) V hiiMiiiy) 

JB ' 4, 



f X7 x G B (x,w) V iiiwMy^ise-'^) 
\V x G B (x,w)\\ J2 U^)\Uv)\e-^ s 



< e^s- 1 



\\V x G B (x,w)\\K B /8 (w,w)?K B /8 (y,y)? 



< e-i^s'i- 1 f \\V x G B (x,w)\\ < e'^^s'l^R, 

JB 



since by (|3.5.6p we have ||VG (x, OIIl^Bj?) -R- If we now take s = t then, by the Euclidean 
estimates and (|3.5.34p . for A\ large enough, we obtain both the lower and upper bounds (13.5.38[) . 

To prove estimate (13.5. 37[) , we use the above estimate and notice that we also have (from Lemma 
13.5.31 and Weyl's Lemma for the ball with Dirichlet boundary conditions) 



<R- d E »iR{hR 2 ) m ° c+1 e- fiiS 

< R-^sR-^-W^R- 1 J2 1 



<A*R-\sR- z ) 

<-(sR 

s 



-2^-2/3,^-1^ 



□ 



Lemma [3.5.101 will be used to get Proposition 13.1.21 for the case of a manifold. We will need to 
improve estimate (|3. 5.371) . which in turn requires the following: 

Lemma 3.5.11. Let \y\ < Rl ) r < =j and < r. Let B v (s r ) be Brownian motion started at y. 
Then 



P( sup \By(s')\ > \y\+r) < d , Cml „. Cmax e 

0<s'<s 



-c (d,c m i n ,c max ) 



Proof. This follows from Lemma f3. 3 .41 □ 

Proof of Proposition \3.L~B: case of Br(z) with metric at least C 2 . By rescaling we may assume that 
R < 1. We upper bound <5o so that \g tk {x) — 5 lk \ < eo where eo is as prescribed by Lemma [3.5.101 
(this is done as in (|3.1.16| ). 

Estimates (|3.1.1[i and the first part of (|3.1.4[) follow from the Euclidean case and estimates 
(I3.5.35|) and (|3.5.36|) . Estimate ([3~Ll]) and estimate (|3~1~3|) follow from (|3.5.38|) and Euclidean ball 
estimates. 
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We now turn to the second and third parts in (|3.1.4p . Without loss of generality we identify 



z = 0. Let a be such that a ^ -k — \- Define stopping times ti,T2, ••• by 

3 = 1 1 

n 

r„=inf{s' :\B z { S ')\=aRY,-}- 

3 = 1 3 

For n > 1, define the set of paths 

B n = {U3 G fi : T n {w) < (1 - 2- n )s} . 

For n > 1 define G„ C B ra as 

G„ = £>„ N -B71-I 

and 

We estimate using Lemma [3.5.111 

P(Gx) <exp(-c'^i? S -5) 

and for n > 1 

P(G„) < exp (^-c'-^-Rs-^ . 

We need another lemma: 

Lemma 3.5.12. The set {u> G l~i : t„ < s V(n. > 1), u> ^ UG n } Zias probability 0. 
Proof. 

(w6Jl:r n <sV(n>l), w ^ UG„} = S O : r„ < s V(n > 1), w ^ U5„} 

= {w6fl:s>T„>(l- 2~")s V(n > 1)} 
C {w £ : t„ - r„_i < 2 _n s V(n > 1)} 

However the set {to £ Jl : t„ — r„_i < 2 _,l s} has probability decaying super-exponentially in n by 
Lemma 13.5.111 □ 

We now continue with the proof of Proposition 13.1.21 case of Br(z) with metric at least C 2 . 
Define H n = G n \ (if^ 1 G l ). We now have a disjoint partition of {lu 6 £1 : r„ < s V(n > 1)} (up to 

measure 0) by the collection {H,}. Set Kf (•, ■) := Iff fl(z) (-, •)■ For |y| > f we have 
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Taking gradient and using equation (|3.5.37|) we get 

oc 



< 



d,c min ,c max ,||g|| 



n=l 
oo 



a J2-(2- n sR- 2 )- 2 ^- 1 s^P(H n ) 



<_ 



< 



.llalU 



,c m i„,c max ,||gj| t 



exp 



n=l 
oo 



n=l 



,a2-V2 



n 00 / .,0-1/2 

< ^(2-"/ 2 <^ J R- 1 )- 4 ^- 2 exp f - c '^±-^2"/ 2 ™ 



1? 



Ilsllo 



where we may replace s with t above, since s < t, and each of the summands is increasing in s as 
long as it is sufficiently small with respect to R 2 (independently of n when n > 1). This proves the 
second and third parts of (|3.1.4[) for Br{z) with metric at least C 2 . □ 

Remark 3.5.13. The proof below makes no assumption on the volume of M, and works for the 
case of M having infinite volume as well. 

Proof of Provosition \3. 1~B for the heat kernel of A4, with metric at least C 2 . As for the Neumann heat 
kernel, the starting point is Proposition 13.3.21 which allows us to localize. We use Proposition 13.1.21 
for the ball P>28 r z (z) with metric at least C 2 . For this proof, we denote by C*2[i?] be the C2 constant 
for the Dirichlet ball case, and set K D (•, •) := K 2 s SaR * (-, •), the heat kernel for the ball B(z,28 R z ) 
with Dirichlet boundary conditions. For s < t, 



\K.(x,y)-K?(x, y )\ = 



+ OO 

5> 



K°_ Tn (x n (oj),y)\T n < s P w (t„(w) < s) 



< 



C 2 [B] 



. m 2 

Ms 

cqn. 1 13.3. 101 



This proves ()3.1.1|) and the first part of (|3.1.4|) (see Remark l3.3.3|) . For the gradient estimates, i.e. 
(|3.1.2p . (|3.1.3p . and the second and third part of (|3.1.4p . 

\\W y K s (x,y) - V V K? {x,y)\\ < ]T ||v,E w \k?_ Tn {x n (u), y)\r n < s] P^ (r„(w) < s) 



n=l 



t 



cqn. 113.3. 101 



_d 5qR z ( °2 " ) 



C^[B],5o,5] 



t"2 



t 



Ma 



giving us C9. By Remark 13.3.31 the exponential term from equation <|3.3.10[) can be made small 
enough so that we obtain estimate (|3.1.2p as well as the second and third parts of (|3.1.4p . □ 
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4 The proof of Theorem [2T2T81 

We remind the reader of Remark 13.1.11 which notes that the proof of Proposition 13.1.21 for the heat 
kernel of Ad, with metric at least C 2 (appearing at the end of section \'S. 5. 3[) . made no assumptions 
on the finiteness of the volume of Ad and the existence of C coun t- 



4.1 The case g 6 C 2 

We appropriately choose heat kernels {K t (z, yi)}i=i,...,d, with t ~ R 2 , that provide a local coordinate 
chart with the properties claimed in the Theorem 12.2.81 

Proof of Theorem \2.2.8\ for g 6 C 2 . Without loss of generality we may assume p = R z = 1, and 
thus, by Remark 13.1.11 we may apply Proposition 13.1.21 Let us consider the Jacobian J(x), for 
x e B c1 r z {z), of the map 

$ := R z d t d/2 (t/R 2 z )<f> . 

By ([3X3)) we have \J l3 (x) - C' 2 (pi, \*~_ V y ]\\ )R^ 1 \ < CgR, 1 . As dictated by Proposition EXU by 
choosing Sq, Si appropriately (and, correspondingly, c\ and cq), we can make the constant Cg smaller 
than any chosen e, for all entries, and for all x at distance no greater than C\R Z from z, where we use 
t = t z = c e R 2 for $. Therefore for c\ small enough compared to C4 we can write R z J(x) = Gd + E(x) 
where Gd is the Gramian matrix (pi,pj) (indepedent of a;!), and |-E"ij(a;)| < e, for x S B Ci r z (z). This 
implies that R z ^ l {a m \ n — Cde)||w|| < ||J(a;)f|| < i?J 1 (cr m ax + Cde)||i>||, with Ca depending linearly on 
d, where <T ma x and cr m i n are the largest and, respectively, smallest eigenvalues of Gd- At this point 
we choose e small enough, so that the above bounds imply that the Jacobian is essentially constant 
in B Ci r z (z), and by integrating along a path from x\ to x 2 in B Cl R r (z), we obtain the Theorem ($ 
and $ differ only by scalar multiplication). We note that e ~ ^ suffices. □ 

We discuss the proof for g E C a , and Ad has possibly infinite volume in Section l4~2l Such proof 
is based on approximation arguments via heat kernels corresponding to smooth metrics on finite 
volume submanifolds. 



4.2 The case g £ C a 

In this section we discuss heat kernel estimates and the heat kernel triangulation Theorem in the 
case when C a . The key ingredient for the proof of Theorem 12 . 2 . 81 for the case of g € C a , are the heat 
kernel estimates similar to those of Proposition 13.1.21 

Before we turn to the proof of Theorem 12.2.81 for the case g € C a , we need one more statement 
about the case g £ C 2 . Consider the following variant of Proposition 13 . 3 . 2l 

Proposition 4.2.1 (Variant of Proposition I3.3.2p . Assume g E C 2 . Let ai £ SJ and R w < 

dist(u>, dtt), or w S Ad and R w < ru(w). Let z and R z be similarly defined. Assume z ^ Br^(w), 
and w ^ Br z (z). For each path B^ (starting at z), we define ti(uj) < T2(w) < . . . as follows. Let 
T\[ui) be the first time that B^ enters B(w, jR w ) (if this does not happen, let ti(cj) = +00). Let 
z\ = B^(ti). By induction, for n > 1 let T n (cu) be the first time after r„_i(w) that B* re-enters 
B{w, jR w ) after having exited B(w, ^R w ), or +00 otherwise. Let z n (u>) = B^(T n ). Lf T n (ui) = +00, 
let T n+ k{uj) = +00 for all k > 0. Then 



K s {z,w) = J2 E * 

n=l 



< s 



Pw(r n < s) 



(4.2.1) 



where 
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Moreover there exists an M = M(c m i n ,c ma x) such that 



Ji-?/; \ / - x _1 / 2 



P(r„ < s) < d ,M, Cmm , Cmax exp{-(n - 1) ( ) (2ms)- 1 - ( ) (2ms)- 1 } . (4.2.2) 



The proof of this Proposition is along the same lines as that of Proposition 13.3.21 

Proof of Theorem 1 2.2.81 for g G C a (M ) with \M\ < oo . Consider a sequence of metrics {gk} QC 2 (M ) , 
with increasing compact supports {M.k\, converging to g in C a (and therefore bounded in C a ), and 
such that gk is uniformly elliptic with constants ^c m j n , 2c max (which is possible since c m i n and c max 
are continuous functions of the components of the metric tensor) . Let K k be the heat kernel asso- 
ciated with gk- Note that for this heat kernel and its gradient we have bounds, from above with 
constants uniform in k for any fixed compact £ away from dM.. We proceed as in the proof of 
Theorem II. 3.1 in |51j . The key ingredients are uniform (in k, for a fixed compact) upper bounds 
on Kk (which follow from Propositions 13.3/21 and l4~.2.ip . and that {Kk} is equicontinuous, which 
follows from the uniform upper bounds on the gradient of Kk (for a fixed compact we have uniform 
lower bounds on R z and R w and estimate (|4.2.2p ). It could also made follow from Stroock's paper 
(Nash-Moser estimates that say the Kk is Holder of order and with constants depending only the 
cllipticity constants)). The proof of Theorem II. 3.1 in 51] then implies that Kk — * Km as k — > +oo, 
uniformly on compacts. Therefore the uniform (in k) bi-Lipschitz bounds on the map 

x -> (K ktt (x, z/i), ...,K k<t (x, y d )) 

on Br(z), imply the same bounds for 

□ 



5 Examples 

5.1 Localized eigenfunctions 

The following example shows that the factors 7i, • • • , 7d in Theorems 12.1.11 and 12.2.11 may in fact be 

d 

required to be as small as Rz ■ 

Let t below be the golden ratio. Consider the domain £7^ as in FigureHJ with Dirichlct boundary 
conditions. We will let z to be the center of the small square. Let and ifij be the eigenvalues and 
eigenfunctions on f^. Fix A > CN . Let = {j : \j < A 2 }. The cardinality of is uniformly 
bounded above by Weyl Lemma. It is also bounded below, since in this case we can easily obtain a 
reverse Weyl Lemma. To see this, let (dfls)A — {x (£ il : d(x,dil) < a/ A}, so that the heat kernel 
estimates in Lemma \3 . 1 . 41 hold for t = bA~ 2 , and observe that 

#{j : < A 2 } > e +b f £ \^(x)\ 2 e~^ 

Jn s \(dn s ) A x z <A 

> e +b ( [ Kf(x,x) f K s t/2 (x,x)e-^) 

\Jn a \(dn s )A Jn s \(on s ) A J 

> e +b (\ - £) 2 (C x r l - C.e-h- 1 ) > e+ b b- l A 2 
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Figure 1 : Example of localization 



where wc choose b so that the last inequality holds, and thus determine a and A, which are chosen 
so that Lemma \3. 1 .41 and Proposition 13.1.21 hold. 

Fix j — 1. The sequence {Xj}s>o is bounded and hence the families {<Pj}s>o and {\7ip S j}s>o arc 
equicontinuous (by Proposition l3.4.lT) , therefore there exists a sequence Sk — > such that ipj k — > </?j, 

V<^ fc — > V(/?j and A* fc — > Aj. We can repeat this argument for any j < liminffe #jj : A^ fc < A 2 } := 
jmax(^4), which is strictly positive and tending to +oo as A — > +oo, by the above. By a diagonal 
argument, we can find a subsequence Si such that for any j < j max (A), <p£ — > tpj, Vy*' — > Vy>j. 
Let us look at some properties of <pj . Clearly, tpj is an eigenfunction for A with Dirichlet boundary 
conditions on f2 - Since (^^i) 2 + (jfn 2 ) 2 is irrational for any ni,n 2 € Z, every ^ is supported 
in either the small square, or the big square. Recall that z is the center of the small square. For 
any j < j max {A) if ipj has support in 5 then ||Vy>j|| > (t/N)~ 2 . Let (5/ be small enough so that 
||V<^' || > (t/N)~ 2 , for all j < j max (A). By choosing A larger than cs(t/2N), where C5 is as in 
Theorem 12.1.11 all possible eigenfunctions that may get chosen in the Theorem will correspond to 
j < jmax{A), and therefore the lower bound for the 7* is sharp. 

See|http : //pmc . polytechnique . f r/pagesperso/dg/recherche/localizatiorue . htm for nice 



demonstrations of the above example. 

5.2 Non-simply connected domain 
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Domain D, point z, closest point to d D, neighborhood to be mapped. 




0.2 0.4 6 08 1 1.2 1.4 0.2 0.4 0.6 0.8 1 1.2 0.2 0.4 0.6 0.8 1 1.2 

Figure 2: Top left: a non-simply connected domain in R 2 , and the point z with its neighborhood to 
be mapped. Top right: the image of the neighborhood under the map. Bottom: Two eigenfunctions 
for mapping. 

6 Appendix 

6.1 Table of Notation/Symbols 

Below is a table of notation/symbols. For each element, we list the first time it is defined mentioned, 
repeat the definition if it is short, or explain what the symbol is typically used for. 
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Symbol Location 

A, A' see ([3~TB| 

b\ see Proposition EHHl 

/?, /3/ oc see Proposition 13.4. II or Lemma [3.5.31 

B r (x), B(x,r) open or closed metric ball of radius r around the point x. 

c see p.l.lOp 

ci, C2, C3, C4, C5, C6 see Theorem 12. 1.11 or 12. 2. II 

Ci, Ca.C^C^ see Proposition [HII2] 

C3, C4 see Lemma 15. 1.41 

C5, 6*6 see Lemma 15. 1.51 

C7, Cg see Lemma l3.1.6l 

Cg see Proposition 13.1.21 or 13.1.21 

^maxi c m i n see ([2.2. ip 

C count see ([2.1. 2p 

<5o, <5i preceding Proposition l3.1.2l 

d the dimension of the domain or manifold 

A Euclidean Laplacian 

A B Euclidean Laplacian with Dirichlct boundary 

conditions on a ball of radius R 
A B Laplace Beltrami operator with Dirichlet boundary 

conditions on a ball of radius R 
Am Laplace Beltrami operator on M. 

71, ...,7d see Theorem 12. 1.11 or 12.2.11 as well as equation p. 1.121) 

S tJ , g lJ Euclidean and non-Euclidean metric tensor 

g r {-,-) see equation (|3.5.1|) 

G r , G r Green function for the Euclidean or non-Euclidean 

ball of radius r (center is some fixed point). In many 
place r is R which is assumed to have value 1. 

/1, £ eigenvalue and eigenfunction of a Euclidean ball 

(usually Br{z) with R = 1) 

fi, £ eigenvalue and eigenfunction of a non-Euclidean ball 

(usually Br(z) with R — 1) 

A, (p eigenvalue and eigenfunction of the whole (Euclidean) 

domain or manifold 

A, <f> eigenvalue and eigenfunction of a the whole 

domain or manifold, with a perturbed metric 

A L , A H A E see pTfi]) . (j3~i~fj)) and (|3.1.10p 

Q, M. a domain or a manifold 

Pi, P2, P3 see Proposition EHHl 

K Dlr ( B \ K B heat kernel with Dirichlet boundary conditions on B, 

context dependent 

K r heat kernel with Dirichlct boundary conditions on a ball of radius r 

(with some fixed center), context dependent 
heat kernel with Neumann boundary conditions, 
context dependent 

K heat kernel, context dependent 

d p K partial with respect to the second variable of a heat kernel 

p see Theorem 12.1.11 and Theorem 12.2.11 

R z written for p, as in Theorem l2.1.1l and Theorem 12.2. II 

rjj see equation (|2.2.2[) 

s a time (for the heat kernel), usually smaller than t (see 

below) 

t a time (for the heat kernel), usually comparable to the 

square of the radius of a ball being discussed 
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